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We consider no-scale inspired supergravity scenarios, where the gravitino mass and related soft 
supersymmetry-breaking parameters are determined dynamically by radiative corrections to an es- 
sentially flat tree-level potential in the supersymmetry breaking hidden sector. We examine the 
theoretical and phenomenological viability of such a mechanism, when including up-to-date calcu- 
lations of the low energy sparticle spectrum and taking into account the latest LHC results and 
other experimental constraints. We (re)emphasize the role of the scale-dependent vacuum energy 
contribution to the effective potential, in obtaining realistic no-scale electroweak minima, examining 
carefully the impact of boundary conditions and of variants of the minimization procedure. We also 
discuss and implement the Bo (soft breaking Higgs mixing parameter) input boundary condition at 
high scale, therefore fixing tan /3(Bq) at low scales. For general high scale boundary conditions with 
Bo, mo, • • • 7^ 0, our analysis provides theoretical correlations among the supersymmetric, soft and 
vacuum energy parameters and related phenomenological consequences at the LHC. For instance, a 
zero vacuum energy at the GUT scale would lead to a decoupled supersymmetric spectrum, together 
with a light standard model-like Higgs boson at the electroweak scale. Given the experimental ex- 
clusion limits, a substantial class of the boundary conditions, and in particular the strict no-scale 
with too = Ao = Bo = 0, are only compatible with a stau being the lightest MSSM particle. Then 
an enlarged allowed parameter space emerges when assuming a gravitino LSP to account for the 
observed dark matter relic density. 
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I. INTRODUCTION 

Supcrsymmctry (SUSY) has imposed itself as the most popular 'beyond the standard model" scenario for many 
good reasons. In addition to appealing extended symmetry principles, it has the potential to solve some of the 
problems raised by the standard model, even though it was not originally introduced for this purpose. It solves 
the hierarchy problem by protecting the scalar sector from unnaturally large radiative corrections, provided that the 
superpartners lie in the TeV range [1-4] . It also predicts the gauge coupling unification [5-8] at a high scale consistent 
with experimental constraints, and provides very plausible particle candidates for the dark matter [9-11]. Last but 
not least, the very structure of the Minimal Supcrsymmctric Standard Model (MSSM) (and of many non-minimal 
extensions) leads generically to the radiative electroweak symmetry breaking (REWSB) mechanism [12-15]. 

The remaining open questions concern mainly the precise mechanism underlying the supcrsymmctry breaking itself. 
Most present viable scenarios assume that a dynamical or spontaneous symmetry breaking occurs in a hidden sector. 
The supcrsymmctry breaking is then transmitted to the visible low energy sector via different mechanisms depending 
on the models. One of the most popular such scenario is when supersymmetry breaking is transmitted essentially via 
the gravitational interaction, in the gravity-mediated models. In an unbroken supergravity (SUGRA) model [16-20], 
the graviton and its supcrpartner, the gravitino, both have a vanishing mass. Once supcrsymmctry is broken, only 
the gravitino gets a mass via the super-Higgs mechanism. Therefore the breaking of local supersymmetry is directly 
linked to the non- vanishing gravitino mass. In a standard SUGRA scenario with a canonical Kahler potential, when 
SUSY breaking is communicated gravitationally to the visible sector the soft parameters are roughly of the same order 
as m 3 / 2 ~ Msusy /Mp, itself expected to be of order the electroweak (EW) scale. In this way one ends up with the 
correct hierarchy between the Planck scale Mp, the SUSY breaking scale Msjjsy and the EW breaking scale Mew, 
although the requirement of the vanishing tree-level potential is somewhat ad hoc. 

In the no-scale models, the basic idea is that the vanishing of the tree-level potential in the hidden sector direction 
can be automatic for an appropriately chosen form of the Kahler potential. Moreover, the value of m 3 / 2 can be 
fixed dynamically by radiative correction stabilization, and is simply related to other soft SUSY-breaking parameters. 
This no-scale approach has emerged quite early [21, 22], and since then it has been regularly claimed to be ruled out 
and resurrected in different forms several times. However, in most of nowadays phenomenological studies, 'no-scale' 
is often a name for just the specific and very restricted boundary conditions on the SUGRA parameters, namely 
mo = Aq = where mo and Aq are respectively the universal GUT scale values of the scalar mass and trilinear soft 
SUSY-breaking parameters, or Bo = mo = -Ao = in the strict no-scale model (Bq being the soft SUSY-breaking 
Higgs mass- mixing parameter). Although well-known to the no-scale model aficionados, it is worth emphasizing here 
that an essential feature of the original no-scale program is the possible dynamical determination of the gravitino mass 
and other related soft SUSY-breaking parameters, that may be realized with the above boundary conditions but also 
possibly with more general SUGRA ones. More precisely the basic framework [21] is to first assume a specific Kahler 
potential such that there exists a fiat direction (moduli) at tree-level, thus also ensuring automatically a (tree-level) 
vanishing cosmological constant. The SUSY-breaking order parameter is the gravitino mass but is not determined at 
the tree-level, i.e. the gravitino mass is 'sliding'. Then the flatness of the potential can be lifted by (non-gravitational) 
radiative corrections originating from the strong and electroweak sectors at the electroweak (EW) symmetry breaking 
scale. These corrections may trigger, under appropriate circumstances, a non-trivial minimum of the potential as a 
function of the gravitino mass, thus fixing the latter. In principle, this picture does not forbid having weak quantum 
corrections to the vacuum energy, as long as those are of order ~ m \/2 [23]- So actually one has rather an 'almost 
flat' moduli direction of the Kahler potential. Overall the mechanism is somewhat similar to the REWSB mechanism, 
but provides an even more direct (and calculable) link between the EW and SUSY-breaking scales. In this way the 
no-scale scenario relates those two scales more dynamically, explaining naturally the hierarchy msuSY Mp. The 
vacuum thus corresponds to a minimum of the potential with respect to the two Higgs fields and a hidden sector field 
z, whose vacuum expectation values (vev) determine respectively the weak scale masses and the gravitino mass m 3 / 2 - 
Thus, the occurrence of a minimum in the z direction is a consequence of the loop-improved effective potential at the 
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EW scale where the tree-level flatness is lifted. For the whole picture to work, one must furthermore assume that there 
are no stronger, purely gravitational (quantum) corrections near the Planck scale Mp, with dangerous vacuum energy 
contributions of the form A 2 TrM 2 , where A ~ 0{Mp) is an appropriate cutoff beyond which quantum gravitational 
effects are non negligible, and M generically the relevant masses of the high scale hidden sector. In fact specific 
superstring models compatible with the no-scale boundary conditions, are known to avoid this problem [23]. On more 
phcnomcnological grounds one may always assume that those issues will ultimately be solved by a fully consistent 
superstring framework, an assumption no more (nor less) problematic than the standard minimal SUGRA picture 
with (assumed) universal soft parameters at high scale. 



On the one hand, most analyses in the past aiming to determine possible no-scale m 3 / 2 minima were conducted with 
definite approximations in the effective potential and sparticle spectrum calculations. Typically, those studies mostly 
used (one- loop) RGE analytical solutions restricted to low tan j3 values, and also typically neglecting the non-dominant 
couplings, non-RG radiative corrections, and other non-dominant terms in the one-loop effective potential, etc, with 
the legitimate aim of determining (semi)-analytical solutions. On the other hand, the situation since those early 
days of no-scale models has drastically changed concerning the elaboration level of (s)particle spectra calculation, 
so that such approximations are hardly considered satisfactory nowadays. There have been of course numerous 
more recent studies of the viability of the mSUGRA subspace defined by the specific no-scale initial conditions with 
more elaborate particle spectra calculations, and updated phenomenological constraints. However, most of those more 
recent studies are generally not addressing the existence of K.// (7713/2) minima. Actually, to the best of our knowledge 
a systematic study of the occurrence of non-trivial no-scale m 3 / 2 minima of a well-defined (RG improved) effective 
potential Vf u u(v u , Vd] rn 3 / 2 ), taking into account its full one-loop radiative corrections, has not been done (though a 
number of special cases, other models, or partial studies of those aspects have been examined in the past or recently 
[24], [25]). In addition, in spite of the above mentioned superstring motivations and different appealing scenarios, in 
the following we will essentially consider a more phenomenological approach. More precisely we shall consider the 
standard mSUGRA parameter models, with more specific boundary conditions like in the strict no-scale models, but 
also more general ones, for which we examine the conditions for the emergence of non-trivial 7713/2 minima at the EW 
scale. This is thus a 'no-scale inspired' but more 'bottom-up' framework, with the goal of determining what kind of 
high scale parameter relations can emerge from our analysis. 



The paper is organized as follows: section 2 gives a short reminder of the basic no-scale supcrgravity scenario. In 
section 3 we specify our procedure for minimizing the loop-improved effective potential, examining some salient features 
to take into account. We emphasize in particular the role played in the V e // (7713/2) minimization by the necessary 
scale-dependent vacuum energy contribution to the effective potential. This point was indeed raised earlier [26-28], 
and in fact our practical procedure is closely related to the latter work. (However at the time those analyses used 
semi-analytical approximations essentially similar to the ones mentioned above, while we will perform more complete 
numerical studies based on the available present SUSY spectrum calculators) . In section 4 we discuss in some details 
important generic properties and results for the minimization of the loop-improved effective potential with respect to 
the extra soft-breaking parameters. We emphasize also the differences in choice of tan /3 or Bq as input parameter, the 
latter being the consistent choice in no-scale scenarios, which has non-trivial technical as well as phenomenological 
consequences. Section 5 examines a few different representative parameter cases, either in the strict no-scale models 
or its generalizations, and following well-defined prescriptions for a loop improved effective potential. The main 
phenomenological constraints affecting the viability of the theoretical results when confronted with collider and other 
experimental limits are illustrated. We also explore the constraints given by the dark matter relic density, either in 
the standard neutralino LSP scenario, or considering alternatively that the LSP is the gravitino, which can be a priori 
assumed in a generalized no-scale scenario. Finally we give some conclusions in section 6. 
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II. BASICS OF NO-SCALE SUPERGRAVITY 



For completeness, we review very sketchily in this section the essential features of the original no-scale models, 
referring for more details to the pioneering literature [21], [29], [22]. The TV = 1 supergravity Lagrangian is fully 
determined by the Kahler potential and superpotential. For simplicity, we focus here only on the matter chiral 
superfield dependence (referring only implicitly to the gauge vector supernelds and gravitino supermultiplet sectors). 
The gauge kinetic function f a b{fa), the Kahler potential K(fa, fa) and the superpotential W(fa) are specified in terms 
of the chiral supernelds fa and their complex conjugate fa. Here fa denotes generically all visible and hidden sector 
supcrfields, possibly transforming under some gauge groups. In terms of the Kahler function 

&{fa,fa) = 5 h In g — (2.1) 

m p m p 

where m p is the reduced Planck scale m p = Mp/y/Sn, one obtains the F-term part of the scalar potential 

v F = mi e G {( ;.(;'(;- - 3) (2.2) 



where 

c - s I' °^£k' GSs(£rv (2 - 3) 

Due to local supersymmetry, the vacuum energy deduced from (2.2) is in general non vanishing even before SUSY 
breaking. After SUSY breaking the gravitino acquires a mass to 3 / 2 given by 

m 3 / 2 = m p e< G >/ 2 (2.4) 

where (G) is a function of the vevs of the scalar components of a sub-class fa; of the chiral supcrfields responsible for 
the SUSY-breaking, whatever the underlying breaking mechanism may be. The corresponding non-vanishing F-term 
vevs F k = ml(e G / 2 G k ) yield the SUSY-breaking mass scale 

Msusy = (FkG^F-) 1 / 4 (2.5) 

Fine-tuning (Vf) to to keep the (tree level) cosmological constant around its observed value fixes (GkG k ) uniquely, 
leading to the well-known relation 

m 3/2 = ^Sfsr . (2.6) 

If SUSY breaking is communicated gravitationally to the visible sector, one expects generically the soft parameters 
m ao ft to satisfy 

m s0 ft = 0(1) x m 3/2 (2.7) 

and m 3 /2 not far from the electroweak scale. Moreover, the generic magnitude (Vf) ~ C( m 3 /2 TO p) ( see (2-2) and (2.4)) 
exacerbates the fine-tuning of the vacuum energy to zero for m 3 / 2 of order the EW scale. No-scalc supergravity [21], 
[30] , [22] was introduced to ensure naturally, through a suitable choice of the Kahler potential, a vanishing potential 
at the tree-level in the hidden sector scalar fields directions at every value of these fields. Since on the one hand the 
0(rriy 2 m p ) magnitude in Vf is now flattened, and on the other the vevs of the hidden sector fields, and thus m 3 / 2 , 
are undetermined at the tree- level, m. 3 / 2 will be fixed at the loop level through (generalized) REWSB which takes 
place in the observable sector, with m 3 / 2 ~ 0(mz) as a natural outcome. This holds only if a large mass scale M 
that can be present in the observable sector, such as a GUT scale, does not contribute to the potential by quantities 
of 0(ml /2 M 2 ). 
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Assuming for simplicity that the hidden sector contains only one chiral supcrficld z (say <fii of the above set of <?Vs), 
the simplest Kahler potential realizing the above potential flatness, entailing a (non-compact) SU(1, 1) symmetry, is 
given by 

K = -3 ln(z + z) . (2.8) 

This Kahler potential has to be supplemented in realistic models by other Kahler potential and superpotential parts 
depending on the visible sector helds. A generalization of the corresponding Kahler function G(z, 4>k) 

G = — ^ + z - + In ( 2 -9) 

where the </>fc's (with k > 2) are all in the visible sector, was found to have very nice properties, either i) when / is an 
arbitrary function but W trivial (e.g. W = m^), or ii) when W is an arbitrary superpotential (e.g. that of the MSSM 
or of an extended GUT model) but / taking the special form f{(f>h,4>k) = I0fc| 2 - In case ii) the SU(1, 1) symmetry 

fc>2 

of (2.8) is extended to SU(n, 1) where n — 1 is the number of fields in the observable sector. A key point for both 
i) and ii) cases is that SUSY-brcaking in the hidden sector leaves the visible chiral supcrficld sector supcrsymmetric. 
This is welcome particularly in case ii) where a full GUT sector can be accommodated, including the MSSM as the 
low energy effective theory, since the non-transmission of SUSY-brcaking will protect the visible sector from the large 
0(m^ 2 MQ UT ) effects mentioned previously. However, for the same reason contributions ©(m-g^m^) will not be 
present either, thus preventing the usual radiative EW symmetry breaking and the ensuing dynamical determination 
of All soft breaking squark, slepton and Higgs masses and couplings are thus vanishing at all scales, and in 

particular the universal parameters at the GUT scale, 



m = A = S = (2.10) 

The only source left for SUSY-brcaking is in the gauge/gaugino sector through the inclusion of a z-field dependent 
non-canonical gauge kinetic functions f a b(z, ...) which are essentially free in a general supergravity framework [16-20]. 
The ensuing soft breaking gaugino mass terms take then the form \m^/2(Gz/ G Z zdf a b/ dz) . One should still assume 
that f a b is chosen such that the heavy GUT gaugino soft masses remain vanishing so that again large unwanted 
contributions 0(rn^, 2 MQ UT ) arc not present. The remaining MSSM gaugino masses are proportional to ITI3/2 as just 
noted. Within the universal gaugino mass assumption we rewrite this relation for later phenomenological use at the 
GUT scale as 

™3/2 = c 3 / 2 m 1 /2 (2.11) 

Equations (2.10, 2.11) define the boundary conditions of the strict no-scale scenario. It is now possible to imagine 
variants to these boundary conditions. For instance supplementing case 1) with a superpotential in the visible sector 
implies a non vanishing Aq. Another variant is to consider Bq ^ keeping the boundary conditions 

m = Aq = (2.12) 

Analogous relations arise naturally in low-energy effective models for some specific string theories [23, 31, 32], 

leading to generalized boundary conditions depending on the type of the string and the compactification mechanism. 

For instance, in the dilaton-dominated SUSY breaking scenario, the conditions become 

1 2 r 

m = —/^ m i/2 > A o = - m i/ 2 , B = -f^ m i/2 1 m i/2 = v3m 3/2 (2.13) 

with thus non-zero values of mo, Bo and other high-scale parameters, but all related to the unique SUSY-breaking 
scale, ~ rrii/2- In contrast, note that the strict no-scale relations in Eq.(2.10) are equivalent to the so-called moduli- 
dominated SUSY breaking superstring model. It is now possible to study the electrowcak potential and to predict 
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values for m\i<i instead of m 3 / 2 by an extra minimization in addition to the ordinary EW minimization driven by 
the REWSB mechanism. Using a simplified analysis with approximations allowing analytical handling of expres- 
sions [21], [22], the preferred values appeared to be of order 

m 1/2 ~ 0(m 3/2 ) ~ 0(m z ) . (2.14) 

As mentioned before, strict no-scale supergravity is characterised by the specific boundary conditions at GUT scale 
given by Eq.(2.10). We are left with only m,\i<z as a free parameter, parameterizing the supersymmetry breaking, and 
driving the other parameters through rcnormalization group evolution (RGE) effects. 

In the following we consider a more phenomenological approach, essentially string model-independent, motivated 
by the fact that the ultimate superstring framework, and even more how it is linked to the GUT scale is not yet fully 
established. We will thus assume the most general standard mSUGRA high scale parameters and boundary conditions, 
but study also the special cases of no-scale (2.12) and strict no-scale models (2.10). This more phenomenological 
approach aims to concentrate more on the conditions for the emergence of a third non-trivial minimum of V e f / with 
respect to m 1 / 2 a t the EW scale, without too strong prejudice on the high scale models. We will assume the following 
generic form of the boundary conditions 

Bo = b Q m 1 / 2 , m = x m 1 /2, A = a m 1/2 ( 2 -!5) 

where bo, ao, xo are mass independent constants taken as input parameters. A further input information about the 
supersymmetric /^-parameter is also needed, even though the value of /x at the electroweak scale will be as usual 
eventually fixed (up to a sign) by the REWSB conditions. Indeed since TO1/2 will be determined dynamically from the 
potential, it is important to know beforehand whether [i has a functional dependence on mi/2, for instance through 
its boundary value ^0 a t the high scale. For reasons which will become clear in the sequel, we will adopt throughout 
the paper the boundary condition 

/xo = c°to 1/2 (2.16) 

with c° an independent constant. Although /j is a supersymmetric parameter, such an assumption is well-motivated as 
there are various mechanisms where it can be related to the SUSY breaking order parameter 771,3/2 within supergravity. 
Typically, this can occur through a non-minimal term in the Kahler potential involving the two Higgs supcrficlds and 
some gauge singlet superfield whose F-term triggers supersymmetry [33, 34], or alternatively through a minimal term 
in the Kahler potential and the addition of an R-symmetry breaking constant in the supcrpotcntial [35] , (see also for 
instance [36] for an early review of other possible mechanisms including the superstring induced ones). 

We close this section by stressing that the usual free parameters of mSUGRA (too, my 2 , Ao, tan/3, sign(/_t)) have 
been traded here for bo, ao, xq and (sign of) c° , thus with one less free parameter, m 1 / 2 , to be determined dynamically, 
hence a more constrained scenario. Furthermore, in the spirit of no-scale, the dimensionlcss parameters bo, ao, xq and 
c° are expected to be of 0(1), or else strictly vanishing. These features are useful criteria distinguishing the no-scale 
scenario from the less constrained mSUGRA, even if both can lead to similar low energy MSSM spectra. 

III. PRESENT SITUATION VS EARLY NO-SCALE MODEL ANALYSIS 

In this section we first review some rather generic and important features of the no-scale mechanism, by which 
a dynamical determination of the soft parameters is provided via the extra minimization of the effective potential. 
Consider the familiar MSSM effective potential, for the moment just at tree-level for simplicity, which reads 

V tree = mj\H u \ 2 +m 2 2 \H d \ 2 - B^H u .H d + ^±il(\H u \ 2 - \H d \ 2 ) 2 (3.1) 

o 

in the relevant electrically neutral Higgs field directions 
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where m\ = m 2 H + fi 2 , m\ = m 2 H + fi 2 , g,g' denote respectively the SU(2)l and U(l)y gauge couplings and the 
'.' denotes the SU{2) scalar product. Once the EW symmetry breaking mechanism occurs, the Higgs fields develop 
non-vanishing vacuum expectation values (h u ) = v u , {hd) = Vd, and the EW extremum is characterized by 

dV dV 

dHjM MH ° )Md={H « } = ' dHdTM MH ° } ' Hd={H « } = (3 - 3) 

where i = 1, ...4 runs over the four field components of the H u and Hd doublets, of which only the two conditions 

dV dV , AS 

oh u dh d 

are not trivially satisfied, and allow to determine Vd and v u in terms of m 2 ,m 2 ,B^, and g 2 + g' 2 . Since the gauge 
invariant point Vd = v u = is also a solution of (3.4), one should require the consistency condition 

m\m\ < (Bfi) 2 (3.5) 

to assure that this point is not a minimum so that the electroweak symmetry is indeed broken. Note that one has to 
require as well 

m? + ml > 2\Bfi\ (3.6) 

to guarantee the tree- level stability of the potential. 1 Of course one should further consider one-loop (and possible 
higher order) corrections to the effective potential and other related radiative corrections to the sparticle masses 
without which the simple tree-level analysis is not sufficiently reliable. Now essentially the main additional feature of 
no-scale models is to seek for an extra non-trivial minimum: 

dV 



dm 3/2 



(3.7) 



or with m 3 /2 replaced by in the above described case of a unique SUSY-breaking scale that can be conveniently 
parameterized in terms of TO1/2, as discussed in the previous section. The parameters in Eq. (3.1) depend non- 
trivially on toi/2 via the high scale boundary conditions, and follow RGE from high scale values e.g. mi{GUT) 
down to mi{EW) values at the EW scale, where the minimization equations (3.4) and (3.7) are required. Before 
examining in more details the extra minimization (3.7), let us first examine some important aspects in defining the 
actual expression for the effective potential to be minimized. As mentioned before, in the early days of no-scale model 
analyses [21], [22] or even a bit more recently [26-28], a number of approximations were used in order to get analytic 
expressions with a rather transparent picture for the behavior of the effective potential and its possible minima. While 
those approximations were legitimate at the time, clearly the situation will change substantially with an up-to-date 
analysis, potentially affecting the existence and location of possible m 1 / 2 minima of the effective potential. We list 
below some of the important features to be taken into account. Rather than giving somewhat blind final results we 
find instructive to disentangle and discuss the different (tree-level, loop level) contributions as much as possible in 
order to pinpoint what contributions are actually responsible for the occurrence of non trivial m 1 / 2 minima. 

• Perhaps the most relevant point concerns the RG invariancc of the effective potential: in principle one would 
expect that the existence of minima is not strongly dependent on the choice of the EW and renormalization 
scales. However, this is a non trivial issue since even the one-loop improved effective potential for the MSSM 
exhibits a rather important scale-dependence in general, due to its intrinsic non-RG invariance unless one 
subtracts a scale-dependent vacuum-cnergy-like term. Although the necessity of including in general such a 
term was established since the work in refs [38-41], we stress that it has a drastic influence on the specific 
mi/2 minimization results, as was indeed pointed out earlier in refs. [26-28]. This will deserve a more detailed 
discussion below. 



There are connections between these two consistency conditions. If (3.5) is violated then (3.6) is necessarily fulfilled. This implies that 
the scale at which EWSB occurs is always higher than the scale at which the potential becomes unstable. Furthermore, when (3.6) 
is satisfied together with (3.4) then the EWSB extremum is guaranteed to be a minimum and (3.5) is automatically satisfied. These 
connections are typical of the MSSM Higgs sector potential and are not valid in a general two Higgs doublet model. They are also 
modified by loop corrections to the effective potential [37]. 
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• In the standard REWSB mechanism, the occurrence of a non trivial EW minimum at some scale Q ew is strongly 
determined by the driving of m 2 H towards its EW scale value, characterized (very roughly) by m 2 H (Qew) < 0. 
Indeed the (one-loop) RGE for mn u reads: 

8 ^ 2 d¥Q m *« = 3Y t 2 ( m Hu + ™% L + ml B + A 2 )- g' 2 Ml - 3g 2 M 2 (3.8) 

where for sufficiently large Yt the first term on the RHS largely dominates. [Note that the trace term 'TrYm 2 ' 
is absent since we assumed mSUGRA boundary conditions at the GUT scale.] One might expect similarly 
that the occurrence of a non-trivial extra minimum of K//(%/2) resembles the REWSB mechanism, therefore 
relying mostly as a first approximation on the running properties and m 1 / 2 dependence of the relevant Higgs 
sector parameters entering (3.1). This is, however, not the case, the detailed mechanism triggering possible 
mx/2 minima being quite more subtle: within the initial conditions (2.15, 2.16) and even when including the 
RGE running of the tree-level potential parameters, (3.1) does not lead to non trivial mx/2 extrema satisfying 
simultaneously (3.4) and (3.7), [except possibly at TOx/2 = 0: but where the extremum is a maximum!]. Adding 
merely a vacuum energy term will already allow for local minima. More generally, as we shall examine 

later on, the occurrence of (phenomenologically relevant) non trivial m.x/2 minima will result from the interplay 
between (3.1), the vacuum energy term and the one-loop corrections to the effective potential. 

• Further influence on the precise location of mx/2 minima comes from the necessary non- logarithmic radiative 
corrections to (s)particle masses, i.e. that are not determined only from RG properties and can indirectly 
affect the effective potential dependence on mx/2- Among those are the field-dependent contributions coming 
from the one-loop part of the effective potential, most conveniently included, at the EW minima, in the form of 
tadpole contributions [42] . Though these are naively reasonably moderate corrections with respect to a tree- level 
analysis, they can have in fact a strong influence on some crucial relations such as those involving /i, m 2 H . and 
the mz mass at the EW minimum. As a result their global effect may shift substantially the mx/2 minima with 
respect to a simple tree-level analysis. In addition, other non-RG radiative corrections are important especially 
for the top and bottom Yukawa couplings and (to a lesser extent) for the gauge couplings. For example in the 
standard procedure where the top (pole) mass is input, one extracts the Yukawa coupling values at some chosen 
input scale Qi n (typically mz or EW) from the relation 

mZp = Y t (Q m )v u (Q m )(l + 6* c (Q m ) + ■■■). (3.9) 

Here the one-loop non-logarithmic (SUSY and SM) corrections S^ c relate the pole mass to the running mass, 
and ellipsis stands for higher order corrections. Now, the supersymmctric and standard model contributions to 
Sy are quite large and positive in most of the mSUGRA parameter space, such that the extracted value of 
Yt{QEw) is substantially smaller than what it would be in a pure RGE approximation neglecting Sy . Then the 
precise value of Y t has an important impact on the mr H running among other things, thus also on the subsequent 
determination of other relevant parameters B, \x via the EW constraints. This emphasizes the importance of 
controlling all sources of radiative corrections for a better determination of the no-scale mx/2(m 3 /2) minima. 

• Finally there are other minor differences between up-to-date standard SUSY spectrum calculations and the 
above mentioned approximations, like the fact that the RGEs are solved numerically for any tan /?, rather than 
analytically for a restricted range of small tan/3 values. In the following we also mainly examine the influence of 
considering consistently the soft breaking Higgs mixing parameter, Bq, to be an input at the high scale, which 
is quite different from considering tan/3 = v u /vd input (at low scale). 
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A. RG invariance and the effective potential 

The tree-level potential (3.1) is known to have in general an unwelcome scale-dependence. A first step to improve 
this situation is to consider the one-loop improved effective potential [43] defined in the DR 1 scheme [44, 45] as 

V eff = V tree + V^ loop = V tree [H u , H d ](Q) + Y^(-l) 2n M*(H u , H d )Qn M ^ H "' H ^ _ |) ( 3 .10) 

n 

where M n are (ficld-dcpcndent) mass eigenvalues and the summation runs over all (s)particle species and possible 
degeneracies due to color, flavor, etc. In the general MSSM the one-loop term in (3.10) takes explicitly the form 2 , 

^i-/ oop = ^M0°) + 2^^^ + ) + 2^^(/)-2 MW)-4X>(C0-16MS) 

0° 0+ / t=l,...4 i=l,2 

-I2{h{t) + h{b)) - Ah{r) + 3h(Z) + 6h(W) (3.11) 



in the notations of [44], where the name of each particle denotes its squared (field-dependent) mass and 

64^ (ln Q2 - 2 



H')^^3^-h- (3-12) 



The expression of RG invariance is formally 



Veff = (3.13) 



where Xi designates generically all relevant couplings or masses, with their corresponding beta functions, and 7^ 
are the anomalous dimensions of the Higgs fields. As is well known the practical cancellations of scale-dependence 
implied by RG invariance generally occurs among terms of different perturbative orders. Thus the (one-loop level) 
cancellation of the scale dependence would be expected to occur between the relevant one-loop beta functions parts in 
Eq. (3.13) acting on the tree-level parameters in Vtree and the explicit Q-dependcncc in the one-loop term above, up 
to higher order (two-loop) remnant terms. However, in general in the presence of massive fields, and in particular in 
the MSSM due to the SUSY-breaking terms, this does not work so, because the effective potential in the form (3.10) 
is not a proper RG-invariant physical quantity, so that there are remnant terms of one-loop order that do not cancel. 
Perhaps rather curiously, apart from early hints [46] this fact was not fully appreciated until the early nineties, where 
different detailed prescriptions were proposed [38-41], all pointing out the necessity to include a 'vacuum energy' piece 
in addition to the above one- loop effective potential (3.10). One simple prescription is to subtract the field- independent 
zero-point (vacuum) energy [41]: 

Vf u ll = Vtree + V\-l oop — Vl-J op,su& (3-14) 

with 

Vi-i 00 p^ su b = Vi-i oop (H u ,Hd = 0) (3.15) 

which can easily be shown to have the same one- loop RG- running as the remnant part from Eq. (3.10), therefore 
cancelling the scale-dependence in the latter up to higher (two-loop) order terms. This subtraction is by construction 
similar to the supertrace in Eq.(3.10) but with a spectrum involving only soft terms and the supersymmetric /i 
parameter 3 . However, this subtraction is only one possible prescription, sufficient for RG invariance properties at 



2 One should in principle also include in (3.11) the one-loop contribution of the gravitino, — 4/1(63/2)- This contribution makes, however, 
very little numerical differences in our analysis, at least as long as m 3 /2 < m i/2> since it has a rather small weight relative to the total 
sum over all contributions to (3.11). For generic input values and = m i/2 it is typically an 0(1%) effect. 

3 For v u = v,} = all particle masses originating from EW symmetry breaking are thus vanishing. The chargino and neutralino masses 
are respectively \M%\, \ft\ and |Mi|. \M%\, \(J,\ with two degenerate ones. Similarly the sfermion and scalar sector eigenmasses depend 
only on the soft terms and the fi parameter. 
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this one- loop order, but having some limitations and unwelcome features [41], such as that the subtracted potential 
may become complex 4 A more general convenient prescription consists in adding a running vacuum energy to the 
potential: 

Vfull = Vtree{Q) + ^l-Zoop(Q) + Kac{Q) (3.16) 

where the running of A vac is determined by requiring Vf u u to satisfy (3.13), which leads to the (one- loop) RGE 
equation, 

Q jQ~ Kac{Q) = E(- 1 ) 2 " M «( i/ - H * = «)■ ( 3 - 17 ) 

n 

A reasonably tractable expression at the two-loop level is also available [44]. The RGE of A vac (Q) has no direct 
influence on other RG parameters, and A vac (Q) behaves at the one-loop level qualitatively exactly like the above 
defined subtraction Vi-i oop . su b, thus canceling the remnant non RG-invariant terms from Vt ree -I- Vi oop in (3.10). 
Moreover, the running A vac (Q) will be uniquely fixed once a choice is made of its boundary value at some arbitrary 
initial scale. In view of later discussions where the boundary conditions (2.15, 2.16) will be assumed, we adopted here 
a notation with the tilde to indicate that the boundary value A° ac can be in principle quite general, including the 
possibility that the dimensionless quantity K^ mc /m\^ 2 be TO1/2 dependent, while a A without a tilde will implicitly 
indicate that the boundary conditions A° ac /m 4 y 2 are chosen 171^2 (TJ3/2) independent. We will refer to the latter 
as the untwiddlcd prescription'. This distinction is important insofar as we are interested in the minima of the 
potential with respect to rrii/2, given that one can always add to A t , ac (<5) any arbitrary Q-independent function of 
TO1/2 without altering the RGE properties while modifying the structure of the minima. Throughout the paper 
we adopt mainly the 'untwiddled prescription' version of Eq. (3.16), but also occasionally illustrate the subtraction 
prescription (3.14), in particular in section IV E. Obviously, the latter prescription is but a special case of the general 
twiddled prescription (3.16), and (3.15) a special solution of the RGE (3.17) with a specific boundary condition at 
some scale Q — Qo- Indeed, taking 

A° ac = —Vi-i ooPtSU b(H Ul Hd = 0, Qo) (3.18) 

as a boundary condition, ensures through (3.17) that 

A„ ac (Q) = — Vi-i oop , su b{H u , Hd = 0, Q) (3.19) 

for all Q. It is also obvious from the form of V\-i oop _ su b that the boundary condition (3.18) is of the twiddled type, i.e. 
^vac/ m i/2 * s a function of m.1/2 when (2.15, 2.16) are assumed. It will prove phenomenologically useful to compare 
the subtraction prescription (3.14, 3.15) with the untwiddled prescription of (3.16). Clearly, these are two different 
prescriptions from the point of view of no-scale, since they differ in the mi/2 dependence of the boundary conditions, 
and thus lead to different Vf u u{mi/2) potentials. 

The vacuum energy being field-independent by definition, has no influence on the EW minimization of the effective 
potential, Eqs. (3.4), so that it can be safely omitted in all related issues. But it can have a definite influence on 
the fate of eventual V{mi/2) minima, contributing non trivially to Eq. (3.7), as we will see in more detail later on. 
In a top-down approach, the running vacuum energy allows to choose different boundary conditions for A vac (Q). 
In particular we shall consider different choices of A vac (QcuT) or A vclc (Qew) an d explore the consequences on the 
existence of m 1 / 2 minima for the full effective potential. This point was already noted and studied earlier [26-28], as 
mentioned in the introduction, but those studies relied on semi-analytical expressions within approximations similar to 



4 This happens in particular in the MSSM case, where for v u = v d = the (would be lightest) neutral and charged Higgs states become 
tachyonic whenever the EWSB conditions (3.5) are satisfied. Note however that these problems may be avoided by subtracting at other 
values of the Higgs fields. 
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the ones mentioned above. Following these authors, we parameterize the vacuum energy contribution at an arbitrary 
scale Q in terms of m\/2 (or equivalently rn^n), in the most general twiddled context, as 

Kac(Q) = n{Q)m{/ 2 (3-20) 

where the running of rj(Q) at one-loop is determined by Eq. (3.17) together with the boundary condition defined e.g. 
at the GUT scale as 7]{Qgut) = Vo- We stress here that, on top of the Qgut dependence, fi(Q) can in general depend 
also on OT1/2 (or equivalently m 3 / 2 )- [we will come back to this point later on when discussing equation (4.5).] 

The improvement in scale (in)dependence of V/ U u (3.16), as compared to Vt r ee + Vi oop alone, is illustrated in Fig. 1 for 
both the untwiddled prescription with mi/2-indcpcndent t]q values, and the subtraction prescription that corresponds 
to an m!/ 2 -dependent initial condition 770- Of course the absolute value of Vf u u depends very much now on the initial 
condition for r\, but this is a constant shift as far as the Q dependence is concerned. Note that, strictly speaking, 
to ensure the RG invariance at one-loop level one should not consider the running of parameters within the one-loop 
expressions in (3.16) [41, 47], since those induce formally two-loop order terms. Indeed, as we have checked explicitly, 
the scale independence of the full effective potential at one-loop is almost perfect when freezing the running of all 
relevant parameters entering the different one-loop contributions, while the formally higher order terms induced from 
those runnings produce a remnant but rather moderate scale dependence visible in Fig. 1. We have checked however 
that such spurious effects remain reasonably small in all cases of our subsequent analysis. In particular they influence 
only very moderately the location of vn.1/2 minima, whenever those exist. 
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FIG. 1. The scale dependence of the effective potential, with or without the vacuum energy contributions, for representative 
input values Bo = .2m 1 / 2 , mo — Ao — 0, mi« = 300GeV and m tov =173 GeV. Dashed line: without vacuum energy 
contribution; full line: Vf u ii from Eq. (3.16) with running vacuum energy and 770 = 10 at the strict one- loop order (i.e. no 
running of parameters inside one-loop expressions); dash-dotted line: same as the previous case but with all running parameters 
at one-loop, showing moderate spurious scale-dependence to be cancelled by two-loop contributions; dotted line: same with 
subtracted vacuum contribution from Eq. (3.15). 



B. The fate of nij/2 minima 



Having a well-defined (one-loop) RG-invariant effective potential, Eq. (3.16), we will now examine in more details 
the behavior of its different contributions with respect to mm. By inspection of the RGE, it is possible to infer that, 
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for generic boundary conditions with a linear dependence upon the soft breaking parameters as given in (2.15), all 
the resulting RG-evolved parameters at the EW scale, in the one-loop RG approximation, will have a similar linear 
behavior, 

Bew = bmi/2, m H u = um i/2' m H d — dm 2 l/2 , A t = a l m 1 / 2 . (3-21) 

The unspecified scale-dependent b, u, d and parameters in Eq. (3.21) have of course complicated expression in terms 
of the original ones in Eq. (2.15), but entirely determined numerically by the RGEs. Note indeed that the linear 
behavior in Eq. (3.21) is obtained even for the strict no-scale boundary conditions in (2.10), because an extra m\/2 
linear dependence is induced from the RGE of the gaugino mass parameters Mi. 

Although the order of the three different minimizations in Eq. (3.4) and (3.7) is in principle irrelevant, before 
examining the mj/j minimization it is much more convenient to sit first at the EW minima, which greatly simplifies 
the procedure. (This is also because we do not consider the full expression of the one-loop effective potential for any 
H u , Hd field values, since it is equivalent and more convenient to put the tadpole contributions when the effective 
potential is evaluated at the EW minimum). This is just the familiar way of expressing the EW minimization Eqs. 
(3.4) as constraints to express Bew and Hew m terms of the other parameters: 

Bfi = (ml a + m 2 Hd + 2 M 2 ) (3.22) 

m 2 Hd - m 2 H J & n 2 f3 g 2 + g' 2 2 
~ tan 2 /3-l V {6 - 26) 

where tan/3 = v u /vd, v 2 = t> 2 + v 2 d in our conventions and rh 2 H . = rn 2 H . + SHi where SH Ut d denotes the corrections 
implied by one-loop tadpoles. The latter have the generic form 

1 ]lf2 

SHi = 16^ c„(-l) 2 "M 2 (ln -f - 1) (3.24) 

n 

where c„ are the different couplings to the respective Higgs fields H u , Hd of all relevant particles in the sum over 
n. Notice that, for reasons that will become clear in section IV C, we did not yet use in Eq. (3.23) the additional 
constraint that the Z mass should be reproduced at the EW minimum, i.e. the extra constraint: 

9 -l±^v 2 ee m| . (3.25) 

As discussed above, one expects the soft-SUSY breaking parameters m 2 H . and B in no-scale scenarios to be directly 
related to the single source of SUSY-breaking, m 3 / 2 (or equivalently mxi'i)- Concerning the supersymmetric fi 
parameter, as stated at the end of section II, it may either be considered as an independent parameter or else related 
to TO1/2 at high scale. But it is in both cases entirely determined, at the EW scale, via the constraints (3.23). 
From (3.22 - 3.24), one obtains, after straightforward algebra, the effective potential at the EW minimum in the form, 

Vf EWmin = _l^J>l v A {l _ 2s 2 } 2 _ ^2 ^ + ^ ^ + + ^ (3 26) 

with cp ee cos/3, sp ee sin/3, and all terms have implicitly a scale dependence here omitted for simplicity of notation. 
Requiring further the constraint of correct physical Z mass (3.25), one obtains 

VfJ mm (m z fixed) = - 2(g2 m J g , 2) (l - 2s 2 ) 2 - 2^L(c 2 5H U + s 2 8H d ) + V loop + A vac (3.27) 

which is formally different from (3.26), in particular as far as the functional dependence is concerned. We will 

come back to this point in more detail in the next section. Note that this difference is strenghthend by the fact 
that the dependence on m# u , m# d has disappeared from the tree- level term of (3.26): indeed, away from the EW 
minimum the (tree- level) potential in Eq. (3.1) depends on the five parameters v u , Vd, m\, m|, B[i (let alone the two 
gauge couplings). Now B and /z can be eliminated upon use of the EW minimum constraints (3.22)-(3.23), but the 
very structure of (3.1) implies that, at the EW minimum, a third parameter disappears so that (3.26) in fact depends 
only on two independent parameters, that we may choose here for convenience to be v u and Vd- 
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IV. LOOKING FOR MINIMA OF THE RG INVARIANT EFFECTIVE POTENTIAL 

In this section, and before entering a more phenomenological discussion, we examine generically the possible exis- 
tence and fate of Vj u u{mi/ 2 ) minima for representative input parameters. We will also specify for this purpose some 
important aspects of our minimization procedure. 

A. Bo input 

We shall first consider an important feature concerning the choice of the input parameters. Most scenarios in 
no-scale models imply a fixed Bq high scale value, in particular in the strict no-scale, Bq = as in Eq. (2.10), or in 
the string inspired case (2.13). However, the by now standard MSSM model-independent procedure is to determine 
Bew together with hew from the EWSB minimization conditions (3.4), not caring usually for high scale values 
of Bq. Even in a more phenomenological framework, it is of interest to perform the minimization rather with Bq 
input, considering tan/3 as dynamically determined rather than an input. Essentially one has to consider Eq. (3.22) 
as determining tan/3 from Bew and the other relevant parameters at the EW scale. It turns out to be a rather 
non-trivial exercise to make such a consistent algorithm. Actually the EW minimization condition (3.4) together with 
tan /3(B) determination turns out to give a fourth order equation for tan/3 (with not always real solutions), but this is 
a rather straightforward part of the derivation. More problematic is that the dynamical tan j3 values thus determined 
at the EW scale is a sensible parameter in all subsequent calculations, and in particular it drastically affects the 
influential top Yukawa coupling (due to the low energy matching relations Eq.(3.9)), which in turn is driving strongly 
the RGE of the B parameter. Therefore, the algorithmically non trivial feature is to get consistent values of Bq and 
tan /3, matching both high and low energy boundary conditions, satisfying the EWSB constraints etc, because of the 
induced effects on the RGE. This has to be solved iteratively and a new algorithm was introduced in SuSpect for this 
purpose. 

We illustrate in Fig. 2 the connection between Bq and tan/3 for different representative values of m 1 / 2 - As one can 
see, rather low tan/3 values are favored, especially for increasingly large Bq. Note that Fig. 2 illustrates only the 
mo = case, but the influence of mo 7^ is not very drastic, only pushing the curves for a given Bq ^ to slightly 
smaller tan /3, in an almost parallel way. The obtained values of tan /3 are essentially located in the range 2 — 25 
(increasing for decreasing Bq and increasing mo) for lOOGeV < m 1 / 2 < O(lTeV). It is still possible to reach larger 
tan/3 values, but for very large m 1 / 2 and small Bq ~ values. 

B. Renormalization scale prescriptions and naive RG-improved effective potential 

Since the full effective potential Eq. (3.16) is approximately scale-independent (up to two-loop order), it is convenient 
to go a step further, following [21], and define a naive 'RG-improved' effective potential by choosing the arbitrary 
scale Qq such that 

V^ 1oop (Qq) = (4.1) 

This choice simplifies largely the analysis since one can then consider only the minimization of Vt re e(Qo) + A-vac(Qo) 
which has a relatively simpler dependence on ra^- The latter combination still embeds one- loop corrections, since 
the dependence Qq(Q) at one-loop is implicit and thus consistently absorbed within Vt re e(Qo) + ^vaciQo)- This 
prescription is partly inspired by the construction of RG-improved effective potential [38-41], where RG-resummation 
properties have been rigorously established for simpler effective potentials with a single mass scale. In the present 
multi-scale MSSM case, a rigorous generalization is not available at present, and choosing simply the scale cancelling 
the loop term is certainly not sufficient to establish correct resummation properties, but at least this prescription is 
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FIG. 2. The Bo — tan/3 connection for various choices of mi/2 (with mo = Aq = 0). 



well-defined and greatly simplifies the analysis. An expression for Qq is obtained from Eqs. (3.10), (3.11) in the form 

n -3/4 PY „ E„(2s« + l)M ? t(Q ) In M^Qo) , . 

Q ° 6 ^ 2 £„(2 S „ + 1)M4(Q ) (42) 

which, at the one-loop level, gives explicitly Qq. [The Qo dependence through M n (Q) on the RHS of Eq. (4.2) gives 
contributions formally of higher (two-loop) order. To include these effects one can still solve (4.2) numerically by 
iterating on Qo-] Nevertheless, the occurrence and location of possible no-scale toi/2 minima are not expected to 
depend too much on this prescription, up to higher order small effects, due to the overall approximate RG-invariance. 
We emphasize that (4.2) is only one convenient choice among many possible prescriptions, and in our numerical 
analysis, as will be more explicit later, we also perform and compare the minimization results when using other 
well-defined renormalization scale prescriptions to quantify the residual scale dependence. In those expressions is 
also implicit the field-dependence of the masses: as mentioned above, minimizing with respect to the H u ,Hd field 
is more conveniently done by canceling the corresponding one- loop tadpoles [42], which is completely equivalent to 
minimizing explicitly with respect to the H Ul Hd fields 5 . All that is performed numerically within SuSpcct [48], 
involving a number of subtleties like necessary iterations etc. 



5 The equivalence only holds when sitting at the EW minimum, as the cancellation of the tadpoles does not give the field dependence 
away from the EW minimum. 
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C. Different m^2 minimization procedures 

Having determined the EW direction from Eq. (3.4), one obtains the effective potential (3.26) in the EW valley, 
VfuU ( m 1/2)1 and then perform the extra minimization Eq. (3.7) along this valley. In this subsection we examine 
in some details different approaches to this minimization procedure, having in mind what can be most easily implc- 
mcntablc within the present tools for the computation of the MSSM spectra. Though we will carry out a numerical 
analysis with SuSpect at the full one- loop level (including dominant two- loop effects as well), it is useful at a first 
stage to examine some approximations to appreciate the main features of the minimization results. 

At first sight it could be very tempting to simply evaluate the full effective potential (3.16) from the output of one's 
preferred SUSY spectrum code, and to perform the extra minimization with respect to m\i 2 numerically. However, 
there is an important subtlety: in the no-scale approach, v Ul va and m\i 2 should be treated as three independent 
dynamical variables with respect to which we seek a minimum of Vf u ii ■ It is only at such a minimum that one requires 
the physical quantities, such as the pole masses of the Z boson, the top and bottom quarks, the r lepton (as well as 
the ones of all other quarks and leptons), to be properly reproduced. Now the point is that in most publicly available 
SUSY spectra codes [48-51] only the v u , Vd two- parameter EW minimization is performed. Accordingly, the mz, nitojf 
and other mass constraints are hard-coded everywhere, which is perfectly consistent in so far as the EW minimization 
conditions (3.22, 3.23) are by definition valid only at the minimum. In contrast, if one is determining possible m 1 / 2 
minima through a numerical analysis of the shape of V^f val (m 1 / 2 ), while requiring beforehand the above-mentioned 
mass constraints, which is the most straightforward way of using the existing codes, then this amounts to requiring the 
correct physical mass values even when not sitting at the physical vacuum. This induces an artificial dependence of the 
soft parameters on m^j distorting the shape of V^ val {m 1 / 2 ) and thus the location or even the possible existence of 
the minima. More specifically, assuming (2.15) for definiteness and fixing mz away from the m 1 / 2 minimum, induces 
through (3.23) and (3.25) a functional dependence of /1 on m 1 / 2 which is incompatible with an (otherwise perfectly 
acceptable) m^-independent boundary condition for fj,, or even with a boundary condition having an m\i 2 functional 
dependence of the type (2.16). Feeding back in (3.26) this fixed-m^ induced ^(mi/2), would lead to an EW valley 
potential Vfff[ val (m 1 / 2 ) that is incompatible with a large class of possible boundary conditions for fj,, and to the false 

9TO1/2 



conclusion that they are forbidden. In practice this would show up in a modification of in the identity 



a j/EWvalf^ „\ _ irEWval , full °^ o% 

^-Jfuu fa/»ri-o^ v w + ^^a^< (4 ' 3) 

that would displace the solutions of dr ^ i Vf^ val { m i/2i A*) = 0. A similar rationale holds for the physical value of 
m^° l p when fixed beforehand. In this case equation (3.9) induces a spurious dependence of Y t (Qi n ) on through 
v u as well as the non-RGE loop corrections 5y C {Qi n )- When fed back in the potential along the EW valley, such a 
dependence will again distort unphysically its shape. 

In order to avoid all these difficulties one can adopt two possible procedures: 

(A) derive an explicit form for (3.7), akin to the explicit forms (3.22, 3.23) derived from (3.3), and solve numerically 
this system of three explicit equations. This permits to use consistently the fixed mass constraints (as they 
stand in the public codes without modification), since one is now sitting at the actual physical vacuum. 

(B) deactivate the physical mass constraints in the codes, determine numerically the minimum of VpJ^ (mi/2) 
along the EW valley, then impose these mass constraints once the minimum is found, to check for consistency. 

Implementation of procedure (B) including the full loop-level effective potential, can be very involved and would 
necessitate a highly non-trivial extension of the standard procedures of the various public codes. Although the EW 
minimization part would not change 6 , the further minimization with respect to m 1 / 2 , (3.7), if performed purely 



the explicit forms at the full one-loop level, (3.22), (3.23), should still be solved numerically due to the tadpoles involving a highly 
non-linear dependence on all parameters. 
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numerically, becomes in general a rather involved task; in the EW valley, one will have to scan a two-parameter space 
( TO i/2iM) ■ requiring the correct fJ,(mz) relationship via (3.23), (3.25), as well as consistency of m v t °^ via (3.9), only 
once the minimum is found. This holds irrespective of whether we assume (2.16) or not. However, as we will discuss 
below, the boundary condition (2.16) will have the benefit that requiring the mass constraints before determining the 
minimum, although in principle wrong, can become a good approximation. This allows a simplified version of (B) as 
follows: 

(C) in this procedure the minimum of V^lf™ l ( m i/2) is determined fully numerically along the EW valley as in (B), 
but the mz and constraints are applied from the outset as in (A). 

Procedure (C) is thus desirable as it provides a simpler picture and is easy to use, not needing drastic modification 
of the present codes. The boundary condition (2.16) has also some practical benefits in the context of procedure (A), 
as it allows to obtain an explicit analytical form for (3.7), which otherwise would not be easily tractable if \x were 
taken as a free parameter. In the sequel we will thus stick to (2.16) in all our study when comparing the outcome 
of procedures (A) and (C). Furthermore, as a consequence of the RGE of fj, which is of the form d/i/dlnQ ~ pt, the 
boundary condition assumption (2.16) leads to 

H(Q) = c M (Q)mi /a (4.4) 

that is valid at any scale Q. Now since fiEW = ^{Q'ew^) * s an y wa y fixed a posteriori by the EW minimization and 
the mz constraint, Eqs. (3.23), (3.25), its high scale value is entirely determined through the RGE from [iew- 
Yet, one has still to verify that fiQ remains of order m\j<i in the spirit of the no-scale framework, as emphasized at 
the end of section II. Indeed, this is generically the case (except possibly for tan/3 very close to 1 which is anyway 
phcnomcnologically largely excluded), given the structure of (3.23) and the fact that [i(Q) does not run much. 

We now discuss in more detail the implementation of procedures (A) and (C) as well as the approximate validity of 
the latter. Using (2.15) and (2.16) and neglecting the running of all masses occurring in the one-loop part of (3.16) 
which is consistent at one-loop level strictly, one can derive from (3.16) an 'explicit' equation for (3.7) to be used in 
procedure (A): 

m U2^—V full ^i/2) = =* V full (m l/2 ) + £(-l) 2 "Af„ 4 (m 1/2 ) + \m\, 2 -^- = (4.5) 

1/2 n 1/2 

where we also used (3.20). Various comments arc in order concerning the validity of this equation. The assumptions 
(3.21, 2.16) are crucial to derive this compact expression, since the first term Vf u ii on the RHS of (4.5) originates 
(upon use of mi/2d mi/2 (m^ 2 ) = 4m^ 2 ) from all terms within Vf u u that scale trivially as n%y 2 , while the supertrace 
term comes from the explicit \n.{m\^ 2 /Q 2 ) dependence within (3.11), see e.g. refs. [26-28]. There is however more to 
it. Since the squared mass eigenvalues M 2 in (3.10) are H u -, i^-ficld dependent, then strictly speaking the scaling 
A/ 2 ~ m 1/2 assumed in deriving (4.5) does not hold in general. There should thus be extra terms in (4.5) due to the 
breakdown of the vn\i% scaling. Such terms are indeed present in general, but they actually vanish when the conditions 
(3.3) are taken into account, i.e. at the true electroweak vacuum. 7 Obviously, such a scaling holds as well at the point 
H u = H d = 0, in particular for the RGE of A l)QC , (3.17). It then immediately follows from (3.20), (3.21) and (4.4) 
that the running of fj does not depend on m 1 / 2 . Thus, the only possible dependence of fj on this parameter would 
come from the boundary condition, whence the last term on the RHS of (4.5). This term, usually not considered 
in the literature, will be important when discussing the alternative subtraction prescription in section IV E. At this 
point equations (3.22), (3.23) and (4.5) form a consistent set of EWSB conditions given the assumptions (2.15), 
(2.16). In a more general context, one would still want to consider models where the (supersymmetric) parameter /zq 



This is not specific to our boundary condition assumptions. For any function V(x> <M with \ an d 4> tw° independent fields, one can 
always artificially treat as being linear in \, <f> = \4>> to recast d x V(x, <f>) in the form of a totaJ derivative V(x, X<M> provided one 
considers only the points where d x V (x, <p) = d^Vix? <t>) = 0- 
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is an independent one at high scale so that the scaling (2.16) does not hold, (reflecting the well-known ^-parameter 
problem). In this case, Eq. (4.5) will have to be modifed in a non-trivial way including additional terms involving 
essentially dVf u u/dfj, that are not straightforward to evaluate for VJ oop . 8 

Equation (4.5), together with (3.22) and (3.23), complete the ingredients of procedure (A). They will have to be 
solved numerically as a non-trivial system of equations giving the values of m\/2, v u and Vd at the minimum of the 
potential. 

We turn now to procedure (C) and discuss the degree of validity of the approximation involved when mz is fixed 
before minimization. To illustrate the case we focus first on the tree-level part of the potential (3.26). Before fixing 
mz, there is evidently a non-trivial m.x/2 dependence in this part, in the EW valley; from equations (3.22) and (3.23), 

4) 



v u and Vd (or v and sV) can be re-expressed in terms of the soft parameters and \i as, 



with 



R 2..2\ 1/2 

1008 201 = 11-4— £-) (4.7) 



i% + 2fi 2 , so that 



2 

V t ree,EWmin = ~ 2 (g 2 + g' 2 ) _ ^ COs2 ^l) ( 4 - 8 ) 

which therefore exhibits an explicit non-trivial m\j 2 dependence via the relations in (3.21). Thus one easily infers 
that Vtree,EWmin ~ m \/ 2 - I R contrast, if one imposes first mz via Eq. (3.25) (as is hard-coded in the public version 
of SuSpect and in other public codes [49-51]), clearly the effective potential has then a different m^-dependence as 
is explicit from (3.27), whose tree-level part only depends on mz and tan/3: it has accordingly practically no more 
dependence on (at least a very mild one as compared to the m^ 2 dependence inferred from the previous analysis 
not fixing mz)- This illustrates the point we stressed previously in this section: fixing mz via Eq. (3.23) beforehand 
induces a spurious functional dependence, /i = /<(mi/2,mz), which overwrites any initial assumption about the fi 
boundary condition and thus modifies unphysically the structure of the minima. This remains of course true beyond 
the tree-level. More generally, it is clear from (4.3) that the reliability of procedure (C) will depend on how far or 
close is the spurious dependence fi(m 1 / 2 ,mz) from the initially assumed model-dependence fJ-{m 1 / 2 ). For instance, 
if fi is taken to be a free miy 2 -indcpendent parameter, then obviously the second term on the RHS of (4.3) should 
be vanishing, while within procedure (C) this is not the case and can even lead to substantial differences for large 
ra 1( / 2 , thus degrading the quality of the approximation in determining the actual minimum. In contrast, the situation 
is much more favorable when the boundary condition (2.16) is assumed. As can be easily seen from (3.23), (3.21) 
and (3.25), the spurious dependence is of the form n(m 1 / 2 ,mz) = ( V-tan^ ™]^ — liig.) 1 / 2 ^ m 1 / 2 in the limit 
m i/2 3> mz, and the functional dependence in (4.4) is properly reproduced at the electroweak scale in this limit. 
Moreover, relatively large m\/ 2 values are favored by the most recent experimental exclusion limits. From this rather 
crude analysis, it is expected that as far as the effect of fixing mz is concerned, the (relative) difference between the 
output of procedures (A) and (C) should decrease like ~ m^/rn 2 ^ for increasing m 1 / 2 . We have performed explicit 
minimizations for rather generic mSUGRA input with Bo, mo, • • • ^ to compare both procedures, and checked that 
the above described qualitative behavior is indeed essentially observed. The corresponding m\j 2 minima that can 
be obtained from the naive (incorrect) procedure can differ substantially for 171^2 ~ mz, while for mi/ 2 > 300 GeV 
the difference between the two procedures decreases to reach about 10% for larger mi/ 2 , as will be illustrated with 



The exploration of this more general case will be pursued elsewhere. 
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concrete examples. The residual 10% difference is not due to the fixing of mz but actually essentially to the effect 
of fixing ?nf°p e through (3.9) which was not taken into account in the above discussion. It is, furthermore, always 
positive in a large part of the parameter space. More precisely, the leading non-RGE one-loop SUSY corrections in 
(3.9) have the generic form, 

where we used (3.21) and (4.4) and where r is a complicated function of 6o, Xq, do of order 1-2 depending on the input 
values. The In dependence translates through the induced m 1 / 2 dependence in Y t into a shift of the form C m i/2 
in (4.5) which in fine accounts for the above mentioned constant relative difference in the values of TO1/2 minima. 

With these features in mind, we will illustrate in most of our subsequent numerical analysis, the results of both 
procedures (A) and (C). 



D. No-scale scenarios and vacuum energy 



We now illustrate minimization results for a representative set of parameter values, adopting the previously described 
prescriptions and minimization procedures, as well as the following three different choices of EW scale: 

1) the 'default' scale, 

Q^ ou,t = (m t - 1 m t - a )V2 (4.10) 
largely adopted in nowadays SUSY spectrum calculations. 

2) the scale Qew suc h that 

Vx-loopiQEw) = (4.11) 

as motivated previously in subsection IV B, so that the expression being minimized is V tree + A vac . This scale 
is to be determined dynamically according to Eq. (4.2). 

3) the scale Qew such that 

Kac(QEw) = (4.12) 

motivated by the requirement of a vanishingly small 'vacuum energy' at the EW scale. 9 Eq. (4.12) has also to 
be solved iteratively, since a different choice of Q ew affects the whole spectrum. 

As it turns out, these three scales are all quite different numerically, so that the comparison of the ensuing results 
is expected to be a reasonably good cross-check of the (necessarily approximate) scale invariance of our minimization 
results. Note moreover that the three choices correspond to dynamically determined scales, being all non trivial 
functions of m,\ji rather than fixed values. 

In Table I we summarize, for mo = A$ = 0,Bq = 0.2mi/2 and 770 = 10, the resulting TO1/2 minima and the 
corresponding values of the Qew scale for these three scale prescriptions, using both the (A) and (C) procedures 
described in section IV C. A first welcome feature is that the existence and values of the minima do not depend much 
on the choice of scale, at one-loop order, as expected if RG-invariance is consistently implemented. We have checked 
that the same rather generic properties are observed for other values of the input parameters. (We comment more on 



Note that this prescription is just a convenient choice. It is by no means intended as a cheap solution to the notorious 'cosmological 
constant problem'. For one thing, at the electroweak scale the true vacuum energy is not given by A, but by the value of Vf u u at 
the minimum, which has various tree-level and loop contributions for non-vanishing Vu,Vd and m 1( i 2 . Then one could rather consider 
a scale prescription such as Vf u a(Q) = 0. But again this is nothing but adjusting the minimum, requiring possibly a proper choice 
of the boundary condition A = ^(Qgut), and certainly not more a solution to the 'cosmological constant problem'. More generally, 
the presently measured (small and positive) cosmological constant is a very large distance observable, and its relation to the very short 
distance vacuum energy computed from a well-defined quantum field theory is another side of the unsolved problem. 
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TABLE I. Values of 7711/2 minima and the corresponding values of QEw{nii/2), f° r m o = Aq = 0, Bo = Q.2m 1 /2, T]o = 10 and 
three different scales, using the two minimization procedures (A) and (C). A conservative intrinsic numerical error of about 
1% is to be added, taking into account uncertainties in the RGE and spectrum calculations. We also indicate for comparison 
the value of Q at the EW border of (3.5). 





Qew = (mti m t 2 ) 1/2 


Vl-loop\Q ew) = 


h- vac {Q ew) = 


EW border 


1) procedure (A): 










Qw(GeV) 


610 


307 


500 


700 


m,i/2{'min)(GeV) 


335 


332 


334 




Vew 


1.1 


-0.6 







2) procedure (C): 










Qw(GeV) 


544 


277 


430 


580 


m\/2{rnin)(GeV) 


297 


299 


300 






0.6 


-1.15 








those checks later in this section.) Another feature is that there is a definite, but rather moderate, difference between 
the results of the two minimization procedures, procedure (A) giving generically slightly higher values of mi/2 than 
procedure (C) (by about 10% here for mi« ~ 300 GeV). This confirms the discussion at the end of the previous 
section where the origin of this difference was traced back to the effect of fixing the physical top-quark mass prior to 
minimization in procedure (C). 

But perhaps the most important feature is that there exists a somewhat narrow window for tj(Qew) where the 
vacuum energy contribution is neither too large nor too small, compensating efficiently the behavior of the tree-level 
and loop contributions such as to produce non-trivial m 1( / 2 minima. For each value of 770 , the precise values of these 
upper and lower critical values of t]ew depend on Bq (mildly) and on mo, and also on the choice of renormalization/EW 
scale. For example the case illustrated in Table I leads approximately to: 

-1-5 < v(Qew) < 1.2 . (4.13) 

with the GUT boundary condition 770 = 10. Specific values of ^{Qew) within this range, corresponding to different 
Qew scales, are given in the table. In fact the lower and upper bounds in (4.13) can be respectively associated to 
the lowest and highest values of Qew consistent with REWSB. For instance, r\ > 1.2 could still lead to a minimum 
of the potential in the m 1 / 2 direction, but would require Qew > 580GcV, which is beyond the EW border (see Table 
I ) so that there is no minimum in the v u , Vd direction. As for the lower bound ~ —1.5, although related to the same 
physical feature, its value is simply dictated by the lower bound Qew > 0(mz) which is a natural requirement. We 
stress that, although corresponding to the same value of 770 at the GUT scale, the various t](Qew) m the range (4.13) 
do not lie on one and the same trajectory of the running rj(Q). The reason is as follows: when the EWSB is fulfilled, 
the values of Y t (QEw) determined from (3.9) with the physical vriF"^ imposed, amount to different and incompatible 
boundary conditions for Y t when varying the Q ew prescription. This leads to a modification of the running of the 
MSSM parameters, in particular rrijj {Q), and so to a modification of the values taken by the beta function in (3.17) 
implying different running rj(Q) trajectories for different Qew prescriptions. 

This point is important to keep in mind when discussing how allowed ranges for rj(QEw) such as (4.13) are mapped 
on allowed ranges of 770 at the GUT scale. We illustrate in Fig 3 the connection between low and high energy 77 values 
as dictated by the RG evolution supplemented with the constraint, but now for a unique EW scale prescription 

Q^e/aw t {m\/2) given by (4.10). In this case the induced difference in the boundary conditions for Y t comes from 
the different values of m 1 / 2 minima. Thus each full- line curve in the figure corresponds simultaneously to a different 
770, a different mi/ 2 minimum and different numerical values for the beta function in (3.17). It is then clear why, in 
contrary to what (3.17) would naively dictate, these curves are not globally shifted with respect to each other and can 
even intersect. Moreover, there is in general no simple one-to-one correspondence between the low energy and high 
energy values of 77 when the physical EWSB constraints are taken into account, since varying the Qew prescription 
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TABLE II. mi/2 minima values for the two minimization procedures and different »7o, with corresponding values of tan f3. Other 
input are Bo = .2m 3/ / 2 ! m o = Aq = 0. The variation of m 1 / 2 minima values corresponds to the various renormalization/EW 
scales, similarly to Table I. We also give some of the phenomenologically most relevant sparticle masses (mq m in designates the 
lightest squark of the first two generations) 



min. procedure 




"71/2 


[min) 


tan /3(Qew) 


c% 


m h 


m x 




1 


mf 


m 


qmin 


m- g 


(A) 


10 


332 


- 335 


6.9 - 7.1 


~ .82 


~ 111 


132 - 


133 


125 - 


126 


692 


- 698 


785 - 


791 


(C) 


10 


297 


- 300 


6.9 - 7.1 


~ .8 


~ 110 


117 - 


118 


112 - 


114 


626 


- 632 


708 - 


715 


(A) 


8 


550 


- 570 


7.7 - 7.9 


~ .86 


~ 115 


227 - 


235 


202 - 


209 


1093 


- 1129 


1250 - 


1293 


(C) 


8 


490 


- 505 


7.7 - 7.9 


~ .85 


~ 114 


200 - 


207 


180 - 


185 


984 


- 1012 


1124 - 


1156 


(A) 


5 


1540 


- 1560 


9.5 - 9.7 


~ .99 


~ 121 


670 - 


679 


552 - 


559 


2789 


- 2822 


3240 - 


3288 


(C) 


5 


1340 


- 1360 


9.5 - 9.7 


~ .97 


~ 121 


579 - 


588 


481 - 


488 


2457 


- 2490 


2854 - 


2894 



would result in a beam of trajectories for each value of 770, which can overlap at low Q. It follows that the lower and 
upper bounds on 770 are somewhat tricky to determine precisely, as they correspond to an envelop deduced from the 
largest allowed range for T}{Qew) at the EW scale. Hereafter we only give a qualitative discussion. An approximate 
range of 770 for which m 1 ^ 2 minima exist, is given by: 

< miQcur) < 15 (4.14) 

and the corresponding largest allowed range at the EW scale is found to be 

-3<#f)<1.7. (4.15) 

Actually, this EW range is obtained from the beam of RG trajectories of tjq = 15 alone, when spanning all possible 
Qew choices. It encompasses all other ranges corresponding to lower 770 values. Fig. 3 is also instructive as regards 
the connection between 770 and m^, given that whatever the precise choice of Qew prescription the latter remains 
of the same order of magnitude as mjw. Accordingly, the location of the mj/j minimum is quite sensitive to rjo, as 
illustrated in Table II for three representative values of 770- In (4.14) the largest 10 < 770 < 15 values correspond to 
TO i/2 J$ 330 GeV approximately. Though it will depend on the precise values of the other input parameters Bq and 
mo, one may anticipate that 770 > 10 is at the verge of being excluded in a substantial part of the other mSUGRA 
parameter space by present collider (LEP,Tevatron and LHC) constraints, as we shall investigate in more details 
in section V. On the lower side of the 770 range, 770 ^ 3 in (4.14) corresponds already to mi {min) > 4TeV, and 
m\ii(rnin) increases very rapidly to extremely large 77i 1 / 2 (3 > ITcV) for lower values of 770- 

This behavior is rather generic and not very strongly dependent on Bq . Varying tjiq has a more substantial effect as 
illustrated on Fig. 3: for a given 770, increasing the ratio mo/mi/2 will result in larger values of possible mi/2 minima. 
Including two- loop RGE effects for the vacuum energy, induces some systematic large shifts of about 20% of the m.1/2 
minima found in Table I, but the results remain qualitatively similar. Pushing it to an extreme, one can extrapolate 
this behavior down to 770 ~ which would lead to mi/2 minima of order 50-80TeV, corresponding respectively to 
large (~ 50) and small (~ 10) tan/3 values. All sparticles have then tens of TeV masses while the lightest Higgs mass 
remains of order 135GeV, a not very exciting scenario for SUSY searches at the LHC. 10 We also observe that the 
sensitivity to the different scale choices tends to increase as 770 decreases (m 1 y 2 increases), in other words the precise 
determination of m 1/ / 2 (mm) somewhat degrades, reflecting the increasing influence of ignored higher order terms. 

In Table II we also give, for each 777,1/2 minimum that we found, the corresponding range of values for tan/3 and 
for some of the most relevant sparticle masses, when different scale prescriptions are applied. We indicate as well 
the corresponding values of c , cf. Eq. (2.16); the latter is entirely determined by the RGE and the EW constraints 



It is amusing to note here that such a configuration, with an obviously severe electroweak fine-tuning problem, comes along with a 
vanishing vacuum energy at the GUT scale. 



21 




1 100 10000 le+06 le+08 le+10 le+12 le+14 le+16 le+18 le+20 

Q [GeV] 



FIG. 3. The one-loop RG evolution of the vacuum energy between the EW and GUT scales, for mo — Ao — 0, Bo = .2mi/2, 
Qew = ( m t 1 m i 2 ) 1 an d three different sets of (t]ew, mi/2) as determined by the minimum of the potential; the corresponding 
rjo values at the GUT scale are explicitly indicated. For rjo — 10, we also show the RG evolution including two-loop effects 
(dashed line) and the mo = .5mi/2 case (dotted line), for comparison. See text for further comments. 



(3.23), (3.25) and comes out of order 1, as expected from the discussion following (4.4). 

An important consequence of the analysis is that the contribution of A vac to the vacuum energy, necessary to define 
an RG-invariant effective potential making it an inseparable part of Vf u u , plays also a crucial role in the occurrence 
of the TO1/2 minima. In the early studies of no-scale minima [21] (where this contribution was not included), the 
parameter choices were mostly such that ^1/2(^3/2) minima were of order mz, so that these minima resulted from 
a fair balance of tree-level and one-loop terms in the effective potential. But given the present phcnomcnological 
constraints on 7711/2 an d the behavior of the tree and loop contributions (including a much heavier top quark mass 
than assumed in the early days) the situation has now notably changed. For instance, using the scale prescription 
(4.10), one finds that for larger m^2, say mi/2 > 300GeV (largely favored by the latest LHC results), the occurrence 
of these minima is essentially driven by the balance between A vac and loop contributions. The latter have a clear 
~ mf, 2 behavior, while the tree- level contributions tend to be relatively suppressed for large m^; see Fig. 4(a) for 
an illustration of this case. However, one should keep in mind that a comparison of the relative tree-level versus 
loop level and/or vacuum energy contributions does not make much sense physically since they are not separately 
RG-invariant. Varying the scale Q shifts parts of Vi-i oop into Vt ree or A„ ac and vice-versa. It is even possible to 
choose the renormalization/EW scale Qew = Qoi^i/T) such that one of these contributions, or some combination 
of them, vanishes identically for arbitrary m 1 / 2 - In the example illustrated in Figure 4(b), all of V\-i oop is absorbed 
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(a) (b) 




m |(2 lGcVJ m^lGeVJ 

FIG. 4. relative contributions of Vtree, A„ ac and Vj oop to Vf u u as functions of mi/2 (for mo — Aq = 0), for two different choices 
of the renormalization/EW scale, using (4.10) in figure (a) , and (4.11) in figure (b). 



in A vac and tree-level contributions modifying consistently their individual shapes in m 1 / 2 - However, once combined, 
they lead to essentially the same value of m\/2 at the minimum irrespective of the scale prescription, as expected 
from the RG-invariance of Vf u u and dVf u u/ dmi/2- Comparing Figs. 4(a) and 4(b) shows indeed a sufficiently good 
numerical stability of Vf u a and of the value of mx/2 at its minimum, against a variation of the renormalization/EW 
scale. 

The dominance of the loop and vacuum contributions with respect to the 'tree-level' ones, as manifest for instance 
in Fig. 4(a) for large m\/2 and the choice (4.10) for the EW scale, should not by itself question the validity of the 
perturbative expansion. For one thing, (4.10) is supposed to partly re-sum the dominant higher order leading logs, and 
is thus in some sense safe (notwithstanding the multi-mass scale difficulties mentioned in section IV B). For another, 
the one-loop and vacuum energy relative contributions are further enhanced due to a large cancellation in the tree-level 
contribution (specially for sufficiently large tan/3), which makes the latter a very shallow function of m 1 / 2 for a large 
class of Qew prescriptions. In that sense it is not so much a problem of uncontrolled higher perturbative orders, but 
a rather accidental very mild dependence of the tree-level on the relevant minimization parameter. 

Finally, the two-loop contributions to the effective potential, though certainly non-negligible in practice, are known 
to remain well under control [44, 45] for the rather moderate values of TO1/2 < 1 TeV that are most phenomenologically 
interesting, and the perturbative validity is not endangered. We have included their dominant 0(a s Yf) contributions 
[52] , [44] and found that they stay at a reasonable level so that the resulting minima are rather stable. For very large 
values of m 1( / 2 well beyond the TeV range, there may be a true perturbativity problem, so that one may not trust too 
much the TO1/2 minimization results. 



E. Alternative subtraction prescription 

We have also analyzed for completeness the structure of the m 1 / 2 minima when adopting the subtraction prescription 
(3.15) rather than the untwiddlcd r/ prescription. In this section we compare the outcome of the two prescriptions for 
a representative set of boundary values for B$, uiq and Aq. For the subtraction prescription, Eq. (4.5) can be recast 
in the more convenient form, 

V fuU (m 1/2 ) i^B" 1 ) 2 " \K(^i/2) - M*(m 1/2 ,v u =v d = 0)}=0. (4.16) 

n 
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Since this prescription provides a special case of the m^-dependent fjo boundary conditions, viz (3.18, 3.19), one 
expects it to lead to a different shape of Vf u a(m 1 / 2 ) than the untwiddled 77 prescription, i.e. not just a constant 
shift, but rather a different structure of the minima in the m 1 / 2 direction. For the same reason the comparison 
between the minimization procedures (A) and (C) of section IV C, carried out in section IV D for the untwiddled 
77 prescription, does not necessarily hold for the subtraction prescription. In particular, as we will see, the very 
existence of TO1/2 minima may now not be necessarily guaranteed simultaneously in both procedures (A) and (C), in 
contrast to what was found in section IV D. We show in Fig. 5 the various dependencies of 77 on 7774/2 taken at the 
electroweak scale with Q fixed to Q d ^^ ult as defined by (4.10) and boundary conditions too = A) = 0, Bq = .2m\/2. 
Following our previously defined conventional notations, 77 represents here either — V su b/m^ , 2 or the (untwiddled) 77. 
In order to understand better the meaning of this plot, we stress that the TOi/2-dcpcndence in 77 has three different 
sources — the boundary condition 770 — the EW scale Qew" C (t71i/2) — and the fixing of 771^ and mF t °^ leading to 
I^Ew(mi/2) and Y t0 p(mi/ 2 ). The two full line curves indicated as 'critical minima' on Fig. 5 are actually the locations 
of the TO!/ 2 minima as determined by procedure (A), that is where Eq. (4.5) is satisfied; the red curve corresponds 
to 77 = — V su b/m^ 2 and the black curve to 77 = 77, taken at the electroweak scale. 11 The outcome of procedure (C) 
is always below those critical curves, due to the shift to lower m 1 / 2 minima by about 10% induced mainly by the 
fixing of 77i(°p 6 , as discussed at the end of section IV C and in section IV D, Tables I, II. It is then easy to trace 
the structure and existence of m 1 / 2 minima by overlaying various curves of — V su b/m^ 2 and 77, corresponding to 
different boundary values of 770 and (2.15), and looking for intersections with the two critical curves. For instance 
the intersections between the 77 critical curve and the set of curves with different Tot ^-independent values of 770 give 
the values of m 1 / 2 and t/ew at the minimum of the potential, as obtained from procedure (A). E.g. 770 = 8,10 
depicted on Fig. 5, lead to m 1 j 2 ~ 550GcV, t}ew — -9 and toi/ 2 — 335GcV, view — 1-1 respectively, see also Table 
II. Varying 770, one can easily read out from the general trend of the monotonically increasing 77 curves as compared 
to the monotonically decreasing critical 77 curve, that there are basically always intersections (i.e. existence of 777^2 
minima from procedure (A)), except for too large or too small values of 770. This reproduces features similar to (4.13) 
and (4.15) albeit here for the specific prescription Q^w^ f° r which the lower critical value of tjew is around 0.7 for 
an mi/ 2 (min) around ITeV, as can be seen from Fig. 5. Moreover, since each 77 curve sweeps out all the region below 
the critical curve, it follows that one will always find a toi/ 2 minimum solution from procedure (C) whenever there 
exists one from procedure (A). This is good news as it guarantees the qualitative equivalence of the two procedures 
even if they differ quantitatively. The situation is drastically different for the subtraction prescription. The reason is 
that the 77 curves of — V su b/my 2 have a shape very similar to that of the corresponding critical curve, rendering the 
intersections rather scarce. One can see this in Fig. 5 where the green (dashed line) curve remains very close to the 
critical curve (red full line) but in fact intersects it at mi/2 ~ 210 GeV which is thus the outcome of procedure (A). 
But after intersecting, the — V su b/m^ 2 does not go enough below the critical curve to match the 10% difference that 
ensures solutions from procedure (C). Indeed, we found that in this case the latter procedure does not yield any m l / 2 
minimum for this specific Qew prescription. In fact, even for procedure (A), the minimum is found in the relatively 
restricted range 300 GeV < Q < 500 GeV. This is not surprising in view of the almost parallel red full-line and green 
dashed curves, which makes the intersection more sensitive to spurious scale dependence from higher order correc- 
tions in this prescription. Other —m~^ 2 V su b curves with different boundary conditions are expected to have the same 
behavior (see e.g. the blue curve in Fig. 5) so that one expects generically that, within the subtraction prescription, 
procedure (C) does not yield solutions when procedure (A) does. Thus, in contrast with the untwiddled prescription 
where a mere 10% effect on the minima is obtained between procedures (A) and (C), the subtraction prescription 
clearly provides an example where the flaws of procedure (C) become severe enough to make it qualitatively unreliable. 

In summary, we learned from this section that: 



Each of these two cases satisfies (4.5), where, as emphasized previously, the latter equation has been derived before imposing the physical 
mass constraints. This is encoded in the fact that drjo / dmi/2 appears in (4.5), rather than dfj/ 'dmi/21 so that this term is non- vanishing 
for fj = —V su i,/m^ 2 but vanishes for fj = 77. This is not in contradiction with the fact that fj has in both cases a non-trivial dependence 
on but only at the eiectroweat scale, as visible on the plot. 
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FIG. 5. the m-L/2 - viQ'sw™ 1 *) correlations. All curves (apart from the blue dotted curve) correspond to mo = Ao — and 
Bo = .2m!/ 2 and Q d ^^ ult = {mi 1 m i 2 ) 1//2 - The t wo 'critical minima' full-line curves give the locations of the m.1/2 minima for 
a given 

v(Q d EW UH ) as determined by procedure (A), Eq. (4.5). (The red curve corresponds to fj = —V su t,/mj, 2 and the black 
curve to f\ = r\.) The black and red dashed-line curves give rj(m,i/2, Q ~ Q^v""'*) f° r different boundary rj , as induced by the 
RGE. Similarly, the green dashed-line curve corresponds to the subtraction prescription. 

1) the subtraction prescription, although theoretically sound, does suffer from an increased sensitivity to higher 
order effects making it in practice perhaps less reliable. 

2) the subtraction prescription might seem more predictive than the untwiddled prescription with the extra r/o 
parameter, but this is an artifact of the specific subtraction at H u = Hd = 0. In general one can subtract at 
other arbitrary values of the Higgs fields with presumably different values of the minimization, thus recovering 
more freedom than naively expected. 

3) although, strictly speaking, only procedure (A) is correct, the comparison between the (A) and (C) procedures 
allowed us to assess better the impact of the various contributions (radiative corrections, physical constraints) 
on the determination of the minima. For instance the effect of fixing mz prior to minimization in procedure 
(C) becomes quickly mild for increasing mi/2, while the radiative corrections when fixing similarly m^° 1 ^ lead 
to moderate but incompressible relative differences of order 10%. 

4) The different choices of scale prescription we considered allowed to check (within numerical uncertainties) the 
expected approximate one- loop scale independence of Vf u u and of the resulting m\/2 minima. Also each one 
of these choices can have its own practical benefit as stated previously, Q^^ ut being the simplest in practice 
since readily implemented in most public codes. 

In the following we will rely exclusively on the untwiddled rj prescription, namely with m^-independent rjo at the 
GUT scale. 

V. COLLIDER AND OTHER PHENOMENOLOGICAL CONSTRAINTS 

In this section we examine the present collider and other phenomenological constraints combined with the theoretical 
constraints from the requirement of non trivial TOx/2 minima. Before going into more details let us start with a first 
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lap regarding the constraints on the relevant high scale parameters Bq (or equivalently tan 0) and 770- Qualitatively, as 
Bo increases (here for fixed Aq = 0), the occurrence of m^ minima is not changing drastically, being mostly sensitive 
to rjo values. But higher Bq implies (for fixed Aq = 0) higher Bew-, an d correspondingly smaller tan/3 (see the 
discussion in subsection IV A and Fig. 2). Therefore the lightest Higgs mass tends to decrease, for virtually the same 
m-y/2 values, with consequently a larger exclusion range in the (toi/2j??o) parameter space. But the present lightest 
Higgs bound can be accommodated even for mo = 0, provided that the miw minima arc sufficiently large as can be 
found for specific t)q values. More precisely, taking rjo < 8 and Bq = 0.2mi/2 one finds approximately nth > 114 GeV 
in agreement with present bounds [53]. Note also that a consistent Bq input tends to favor 5 < tan/3 < 25, as long as 
TO i/2 ^ 1 TeV. We give for illustration in Fig. 6 the full sparticle spectrum obtained for the case r/o = 8, Bq = .2m.i/2 
(and mo = ^0 — 0) with a minimum at about mi/2 ~ 500 GcV. The lightest Higgs mass is m/j ~ 114.3 GeV and this 
spectrum passes all other present constraints (including b — > sj, — as well as the recent LHC constraints [54], [55] 
), except for the important fact that the f is the LSP. This is the case more generally in a large part of the parameter 
space. We shall discuss in the next section how to evade this dark matter issue by assuming that the gravitino is the 
true LSP and the stau the NLSP. The longstanding cosmological issues related to the gravitino including gravitino 
LSP as a possible dark matter candidate [56], were also considered in the no-scale framework [29, 57]. 



1000 



H, H+, A 



100 L 



"0,3.4 



"+,1 Xn 



n 0= A 0= B 0= 2 m i/2 

T| =8 m ]/2 (min)=500 



FIG. 6. A typical representative spectrum with a consistent mi/2 minimum. 



A. The constrained MSSM and the LSP issue 

In standard mSUGRA, substantial parts of the parameter space where the lightest ncutralino is the LSP and 
assumed to be the dark matter, can be excluded by the previously established collider and relic density combined 
constraints [58], [59]. In such a scenario, most of these results may be roughly applied in our case, provided that 
one superimposes on those constraints the specific (m^, tan/3) values theoretically constrained by the no-scale Bq 
input, together with the extra constraints on the vacuum energy via 770- This nevertheless deserves a specific study 
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and update, as we illustrate for a few representative cases below. Now in fact, for a large part of the relevant no-scale 
parameter space, the LSP is the charged stau, which is in conflict with the requirement of an electrically neutral dark 
matter. The neutralino mass grows faster with than the stau, which consequently tends to be lighter. Having 

a non vanishing too raises the initial value of the scalars and delays the moment where radiative contributions raise 
the neutralino mass, giving more easily a neutralino LSP. For rather small values of TO1/2, say less than about 300 
GeV, and rather low values of tan/3, there is a window in which the neutralino and stau neutrino are lighter than the 
stau, but this part of the parameter space is largely excluded at present by other lower sparticle mass bounds from 
the LEP [53], Tcvatron[60], and the latest LHC limits [54], [55] . In particular both the right selectron, the lightest 
chargino, and the gluino recently constrained by the LHC, can easily be too light. 

This issue can be solved by assuming that the gravitino is the true LSP, thus lighter than the stau, so that the 
latter is decaying to a gravitino plus a r lepton in the early universe. In our generalized no-scale framework the 
gravitino mass is required to be 0(7711/2) but otherwise essentially free. For this alternative scenario new and specific 
constraints arise mainly from relic density confronted with WMAP results, including also in this case the very relevant 
gravitino thermal contributions to the relic density. 

We thus assume that the gauginos and the gravitino are linked in a manner that allows a light gravitino, though 
not exhibiting a precise relation. Note that the gravitino does not need to be very light, but just slightly below 
the stau and neutralino masses still with J7i3/2/%/2 ~ On purely phenomenological ground we assume the 

simple relation (2.11) and explore the different constraints for representative choices 0.1 < c 3 / 2 < 1. More elaborate 
model-dependent relations are also possible [22], in particular in superstring derived models. 

From the point of view of dark matter relic density constraints, assuming a gravitino LSP opens up a whole new area 
of parameters that is otherwise excluded for a neutralino LSP. Requiring in addition non-trivial no-scale mi/2 minima 
and consistent Bq input, gives tighter constraints on TO1/2 and also on rjQ compatible values as we have illustrated 
above. 

In addition other more indirect phenomenological constraints on supcrsymmctric models, such as those obtained 
from the muon anomalous moment and the B decay observables, will be taken into account. But indirect constraints 
arc less drastic in general, since they can always be fulfilled by additional contributions or slightly modified scenarios. 
While relic density constraints are the most important from our perspective (on top of direct collider limits), since it 
is the only way to put constraints on the gravitino mass in this general no-scale scenario. 

1. Sparticle mass limits 

We use present collider limits on sparticle masses some of which arc, however, model-dependent. For instance 
many available bounds assume a neutralino LSP, or when assuming the decay to a gravitino, that m 3 / 2 is very small 
(^3/2 ^5 1 kcV). It is thus difficult to read out general bounds on the masses from the existing limits. For simplicity 
we may apply conservatively the pre-LHC limits for the colorless MSSM particles that are strictly speaking valid only 
for a mSUGRA scenario with neutralino LSP [53]: 

• neutralino mass: to x o > 46 GeV; this bound could be easily evaded, however, by relaxing the gaugino unification 
assumption, that we took for simplicity but which is by no means mandatory in a no-scale framework. Indeed 
there are little theoretical constraints on the form of the gauge kinetic function which generates a non- zero m.1/2 
value, so that it may not need to be universal [22]. 

• chargino mass : to ± > 104 GeV for mc, > 200GcV; this lower bound can be somehow evaded in the case of a 
very light sneutrino (which contributes to the t-channel with destructive interference). 

• stau : TOf > 86 — 95 GeV, valid as long as (\mf — m x \ > 7 GeV) 

• smuon : to^ > 95 GeV 
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• squarks and gauginos: 

Recently CMS [54] and ATLAS [55] have put new exclusion limits on squark and gluino masses from the LHC 
run at 7 GeV center-of-mass energy with integrated luminosity of 35pb _1 . In the mSUGRA model this translates 
into lower bounds in the m 1/ / 2 _ '7io plane that are more severe from searches for multi-jet events than for events 
containing two or more leptons in the final state. Furthermore, these constraints are not very sensitive to Aq 
and tan /3 as far as the latter is not very large (a regime not relevant to our case) . These limits lead to a lower 
m 1 / 2 bound of about 300 GeV for moderate m values 0-200 GeV. While finalizing this paper a very recent 
ATLAS note [61] appeared, extending the study of jet events with missing energy for a higher luminosity of 
165pb _1 . This appears to exclude < 450GeV for low too values and Ao = tan/3 = 10, thus ruling out 

a priori our first benchmark study in Table I, if applying conservatively the limits valid for a neutralino LSP. 
In any case this is not much of a problem in so far as our different benchmarks are only theoretical examples 
for the occurrence of no-scale minima. Note that the spectrum e.g. shown in Fig. 6 is still (slightly) above the 
border of present exclusions. 

However, it should be emphasized that all those constraints do not necessarily apply in a large part of the parameter 
space considered here, where to^ < m x o , with a gravitino LSP and f i NLSP. Indeed, for the range of gravitino 
masses that we will consider, the stau NLSP is sufficiently long-lived not to produce a signal with missing energy in 
the detector, so that limits obtained from searches for such signals do not apply. Very recent results by ATLAS [62] 
dedicated to searches for long charged tracks give bounds on stable staus NLSP, mf > 136GcV, for a representative 
gauge-mediated SUSY-breaking (GMSB) model. Since this limit concerns a different model, its implication on no-scale 
would necessitate a detailed study. Nevertheless, naively fitting this f mass bound for too = Aq = and an appropriate 
value of tan/3(I?o), leads to a quite similar MSSM spectrum corresponding approximately to mwj > 360GeV. 

2. Higgs boson mass 

Strong limits on the lightest Higgs mass are obtained from LEP and the Tevatron [53]. In fact, apart from a very 
small window for relatively small TO1/2, which is now essentially excluded by Tevatron and the above mentioned very 
recent LHC limits, in most of the no-scale parameter space we are in the Higgs-decoupling scenario tua S> such 
that the standard-model like limit on rrih essentially applies. Thus in most of the parameter space we are considering 
here the LEP II constraint [53, 63] should hold : 

to/j > 114.4 GeV (5.1) 

where the limit depends on tan/3 to some extent. This implies rather stringent lower bounds on m^2 prior to LHC 
squark/gluino mass limits, specially for low tan/3 (i.e. large Bq). Even for Bq = 0, requiring Eq. (5.1) corresponds 
to m 1 / 2 > 400GeV e.g. for ?n = A = 0. Otherwise, the most conservative limit is rrih > 92.8 GeV [53]. To 
take into account theoretical uncertainties, we will use a conservative limit allowing for about 3-4 GeV of theoretical 
uncertainties as is customary. 

3. Muon anomalous moment 

Supcrsymmetric particles can contribute at the loop level to the muon anomalous moment — (g — 2) M /2 [64-69]. 
It could explain the deviation measured [70] : 

Aa M = a e * p - = (22 ± 10 to 26 ± 9) x 10~ 10 (5.2) 

which is about 2 to 3 standard deviations from the theoretical prediction of the standard model (for a review see [71]). 

The contributions to a M in no-scale models is of course a particular case of mSUGRA general contributions. These 
can be quite important when slepton masses are light. Corrections in a general MSSM come dominantly from loops 
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with a chargino and a muon sneutrino and loops with a neutralino and a smuon [66]. The MSSM correction is 
proportional to tan/3 and its sign follows the sign of fi, thus favoring a positive sign for the /i parameter in view of 
(5.2). For not too large tan/3 values and choosing /j, > 0, one can accomodate the preferred range (5.2) for some 
regions of the mSUGRA parameter values compatible with no-scale models (provided of course that to^/2 is not too 
large), as we shall illustrate in the next section. 

4- b — > s 7 

Another largely studied probe for supersymmetry is B-meson physics, in particular the decay b — > which has been 
extensively measured with good accuracy and is theoretically well under control. Indeed, theoretical calculations from 
standard model contributions have been now performed at the next-to-next-to- leading logarithm (NNLL) order [72], 
including also nonperturbative corrections [73]. Confronted with recent experimental measurements [53, 74]: 

Br(b s 7 ) = (3.55 ± 0.24 ± 0.09) x 10" 4 (5.3) 

it results in a discrepancy with the standard model slightly above one standard deviation, therefore potentially very 
constraining for new physics. The possible contributions from the MSSM are dominated by one-loop effects from 
chargino plus stops, and top plus charged Higgses. NLL SUSY-QCD corrections have also been calculated [75]. In 
mSUGRA, contributions can become sizeable for relatively large tan/3 and sufficiently small mi/2, which is not much 
favored in no-scale scenarios, due to the m 1 / 2 — Bew — tan /3 correlations inducing rather moderate tan /? < 10 values 
for not too large mi/ 2 (see the discussion in subsection IV A). We therefore anticipate that b — s- s-f constraints are 
relatively marginal for a large part of our no-scale inspired scenarios, specially for the strict no-scale with mo = 0. 

In practice we have used in our analysis the bounds (5.3) conservatively augmented by theoretical uncertainties as 
quoted e.g. in [72, 73]. The precise limits are, however, not very crucial for our analysis, since as we will illustrate 
they give anyway (mild) constraints only for rather small TO4/2 < 250 — 300 GcV which, for the range of tan/? values 
under consideration, are largely superseded both by the lightest Higgs mass limit from LEP and by the Tevatron and 
the recent LHC limits on mi/ 2 - (In fact for Aq = and too 7^ moderately small, the lightest Higgs mass bound 
irih > 114.4 GeV generally supersedes the b — > constraints even for large tan /3 values up to tan /3 ~ 40 — 50, anyway 
unreachable in the no-scale framework, b — > sj can be much more constraining if — Aq is large enough, |j4o| ~ 1 — 2TeV 
such that the stops could be light enough, see e.g the discussion in [59]). 

5. Dark matter relic density 

The lightest neutralino as a candidate for dark matter has been extensively studied in many scenarios [58], [59]. 
The part of mSUGRA parameter space giving a stau LSP should be normally excluded. This has been one argument 
advocating against the viability of the strict no-scale scenarios. But if one considers a gravitino lighter than the stau 
then this part of the parameter space regains interest. The gravitino dark matter candidate has also been quite studied 
in the past decade. In the case of gauge mediation supersymmetry breaking scenarios, this particle is naturally the 
LSP and was considered both for cosmological issues and in colliders signatures [76-82]. The gravitino can also be 
the LSP and a very interesting dark matter candidate in the context of mSUGRA scenarios [82-95]. 

In the present analysis we will illustrate a few scenarios for the most representative no-scale cases as studied above. 
A more complete study of the constraints obtained for the full mSUGRA parameter space will be done in a forthcoming 
analysis, where we also consider in detail some implications and constraints from Big Bang Nucleosynthesis on such 
LSP gravitino in no-scale scenario. 

For the relic density, we use micrOMEGAs 2.0 [96], to compute the relic density of the neutralino or stau MSSM 
LSP. For scenarios with a gravitino being the real LSP (neutralino or stau being the NLSP), all supersymmetric 
particles decay to the NLSP well before the latter has decayed to a gravitino, because all interactions to the gravitino 
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are suppressed by the Planck mass. We first compute the relic density ^nlsp^ 2 the NLSP would have if it did 
not decay to the gravitino. Then assuming that each NLSP with mass tonlsp decays to one gravitino, leads to the 
non-thermal contribution to the gravitino relic density 

fiNTP^ = Z_V__n NLSp/l 2 (5 .4) 

1 Tl NLSP 

with h = 0.73t!j;o3 the Hubble constant in units of 100 km Mpc" 1 s _1 . 

The gravitino can also be produced in scattering processes during reheating after inflation [57, 97-100]. Following 
[98, 99], the resulting gravitino yield from thermal production is controlled by the reheat temperature Tr as follows, 

« no - ± M (. + >» ) (^) (») 

where i sums over gauge groups, y.;/10~ 12 = (0.653,1.604,4.276), h = (1.266,1.312,1.271) and gauge couplings and 
gaugino masses are calculated using the one-loop RGE. Assuming standard thermal history without release of entropy, 
the gravitino relic density from thermal production is 

njf 2 h 2 = m 3/2 Y 3 T ^(T Q ) S (T )h 2 1 p c (5.6) 

with p c /[s(To)h 2 ] = 3.6 x 10~ 9 GeV and To the background temperature. 
Comparing the total gravitino relic density 

n 3/2 h 2 = njf 2 h 2 + n^h 2 (5.7) 

to the one inferred from the measurements of the CMB anisotropics will constrain Tr and the MSSM parameter 
space. The 3-year WMAP satellite survey has given at 3cr confidence level [101] 

n% M h 2 = 0.105±8;{gS (5.8) 



B. Combined Constraints 



We can now combine both the theoretical no-scale constraints (i.e., the existence of non-trivial mi/ 2 minima and 
the So — tan /3 relationship) and phenomenological constraints on m.1/2 and other parameters from direct and indirect 
search limits at colliders and from other obscrvablcs. LEP, Tevatron, and latest LHC constraints tell us that rather low 
rrc.x/2 are now essentially excluded, also to avoid the 114 GeV lightest Higgs limit and other indirect phenomenological 
and theoretical exclusions. WMAP will also give stringent constraints excluding essentially all low mo values (in 
particular the pure no-scale mo = case) if the neutralino is the true LSP, because the corresponding relic density 
comes out below the observational bound (5.8). 



1. Collider and other phenomenological constraints 



In Fig. 7 we give the present constraints from direct sparticle search limits and the low energy constraints from 
(5.3), (5.2), in the (m^, tan/3) plane for mo = A$ = 0, most relevant to the true no-scale scenario. We do not 
put explicitly the above mentioned latest LHC exclusions [54, 55] on this and other subsequent plots, since these are 
anyway debatable given that in most of our scenarios the neutralino is not the LSP as discussed previously. A number 
of conclusions may be easily drawn from this figure: 

• Due to the Bq no-scale input, one obtains for each Bq value specific mm, tan/3 correlations. In particular for the 
strict no-scalc (2.10) the values of tan/3 are restricted to be > 20 when taking into account other constraints. 
This can be consistent with the Higgs mass lower bound of ~ 114 GeV and falls into the preferred b — > sj range, 
but in all this region of parameter space the stau is the MSSM LSP as indicated, so that a gravitino true LSP 
becomes a very appealing scenario. 
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• For larger Bq values, there can be regions where the neutralino is again the LSP, typically for Bq > .15mi/2 
and sufficiently small m 1 / 2 , see the figure. But this is generally not compatible with the light Higgs mass 
limit, even when allowing a large theoretical uncertainty. In principle there could be a tiny region for such 
Bo ~ 0.2 — 0.5m!/2 values, where for sufficiently small m^ 2 one is no longer in the decoupling limit, i.e such 
that ttia is light enough and the bound in (5.1) no longer applies. However in that case the very recent direct 
limits from the LHC [54], [55] exclude virtually all of this small corner. 

• The requirement of non-trivial mi« minima leads to constraints e.g. in the plane (Bq^tjq) or equivalcntly 
(tan /?, rjo) for mo = 0, or more general ones for mo ^ 0, that we do not give explicitly. Suffice it to say that for 
any 770 in the range (4.14) one can find TO1/2 minima, but present lower bounds on TOx/2 exclude accordingly 
Vo ^ 8-10, approximately. The most phenomenologically interesting range, obtained for not too large TO1/2 < 1 
TeV, corresponds to 5 < 770 8, while for smaller 770 the corresponding TOj/ 2 minima increase very fast as 
discussed before. 




FIG. 7. Direct collider and other indirect constraints in the (rrei/2, tan/3) plane for mo = Ao — 0. The lines for three different 
representative Bo input give the m 1 / 2 ~ tan/3 no-scale correlations. Dark and light grey colors indicate respectively the areas 
excluded by inconsistent EWSB ('tachyonic'), and by direct sparticle mass limits from LEP and Tevatron. The light blue zone 
is excluded by b — > S7 constraints (5.3), while the yellow zone corresponds to values falling into the measured (? M — 2 range of 
(5.2). The dashed lines give the lightest Higgs mass limits and the dot-dashed line the mf = border. 



2. Relic density constraints and gravitino LSP 

In Fig. 8 we show the relic density values in the (mi/2, tan/3) plane corresponding to Fig. 7, calculated for f MSSM 
LSP after it has decoupled from the thermal bath. We also show, in yellow, the small region where Ni becomes the 
MSSM LSP. In the latter region the Ni abundance remains too small to be consistent with WMAP, thus excluding Ni 
as a dark matter candidate but still allowing it as an NLSP with, in this case, a rather small non-thermal contribution 
to the gravitino relic density. The relic density values obtained in the (largely dominant) region where the f is the 
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MSSM LSP only make sense if the gravitino is the true LSP, and will lead in some parts of the parameter space 
to a substantial non-thermal contribution (5.4) to the total gravitino relic density. We illustrate the gravitino total 




FIG. 8. r relic density values at decoupling, in the (mi/2,tan/3) plane with mo = Aq = 0. In all blue regions f is the MSSM 
LSP. Different levels of blue correspond to different values of the f abundance, as indicated; in the yellow region Ni becomes 
ligher than f but with very small relic abundance (see text for more comments), other captions as in Fig. 7. 



relic density, obtained from (5.4 -5.7), in Fig. 9 for m 3/ / 2 = 0.1mi/2 and in Fig. 10 for a higher m 3 / 2 /wi/2 ratio, for 
different values of the reheating temperature Tr. The first value is such that the gravitino is the true LSP in most of 
the parameter space, while in the second case the uncacceptable region where m 3 / 2 > 771 f is enlarged for large tan/3. 
Now one can see that it is easy to recover consistency with the WMAP relic density constraint in a large part of the 
parameter space, provided that Tr is sufficiently large, Tr > 10 6 GeV. There is not much qualitative differences for 
the two illustrated m 3 / 2 masses. In fact we observe that a phcnomenologically most interesting case for potential 
early discovery at the LHC, namely for not too large miw > 400 — 500, and consistency with the — 2, b — > sj and 
WMAP constraints, implies a relatively large Tr > 10 8 -10 9 GeV, the darker red region on the figure. A large Tr is 
also welcome by other independent issues such as thermal leptogenesis scenarios [102]. Indeed, a comparison of Fig 7 
and Fig. 9 shows the interesting fact that a part of the strict no-scale model Bq = tjiq = Aq = is not excluded: there 
is a range, for m 1 / 2 ~ 400 — 800 GeV, tan/3 ~ 20 — 25 compatible with rrih > 114 GeV, the — 2 deviation (Fig. 7) 
and other constraints, provided the reheating temperature is 10 8 — 10 9 GeV. 



3. Generalized no-scale scenarios 

In this subsection we consider one typical example illustrating more general cases with Ao,mo ^ 0. Since Aq has 
little effect on the existence of non-trivial no-scale minima, we illustrate for simplicity here only the case Aq = 0. In 
Fig 11 we show the results of the same analysis as before but for mo = .5mi/2. An important difference with the 
mo = case in Fig 7 is a larger region excluded by b — > sj, which includes a substantial part of the Bq = line, for 
rather small m\/2- Now coming to the relic density constraints on Fig 12 there exists an interesting region (in green) 
where the neutralino, if it is the true LSP, is compatible with WMAP. Note however that this region is almost entirely 
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FIG. 9. Gravitino DM relic density values consistent with the WMAP constraints (5.8) in the (mi/2, tan/3) plane, assuming 
m 3/2 = O.lmi/2 and mo — Aq = 0. The different levels of red correspond to different reheat temperature values. The small 
region where the gravitino is not the LSP is also indicated; other captions as in Fig. 7. 




FIG. 10. Same as Fig 9 but for m 3 / 2 = §wii/2- 



excluded by b — s- for the Bq = line as can be seen on Fig. 11 and also basically excluded by the latest LHC 
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limits on m\/2- But one can easily find appropriate too and -Bo values such that the b — > sj and LHC constraints 
are compatible with a neutralino LSP. A full scan of the parameter space will be explored elsewhere. Alternatively, 
the gravitino LSP case with its relic density is illustrated in Fig 13, where the main difference with the mo = pure 
no-scale case is that the consistent region with il^^h 2 < VlnMh 2 is shrinked to much smaller m 1 / 2 values for which 
there is accordingly a tension with the latest LHC m\/2 lower limits. For sufficiently small (but not yet all excluded) 
TO4/2, one can have the right relic density with a high reheating temperature almost independently of tan/3 similarly 
to mo = 0. 




FIG. 11. Direct collider and other indirect constraints in the (mi/2, tan/3) plane for mo — .5m 1 /2, Ao — 0. See Fig. 7 for 
captions. 



VI. CONCLUSIONS 

We have re-examined generalized no-scale supergravity-inspired scenarios, in which the gravitino mass and all other 
connected soft supersymmetry-breaking parameters can be dynamically determined through radiative corrections, 
triggering a non-trivial minimum of the RG-improved potential. For representative high scale boundary conditions on 
the minimal supergravity model, we have examined critically the theoretical and phenomenological viability of such 
a mechanism in view of up-to-date calculations of the low energy supersymmctric spectrum, taking into account all 
important one-loop radiative corrections. We also have investigated the impact of different prescriptions and possible 
variants of the minimization procedure, paying attention to the extra TO^-dependcncc, genuine or fake, induced by 
the implementation of physical mass constraints with various sources of radiative corrections. 

We emphasize the importance of using a RG-invariant effective potential including consistcnly a scale dependent 
vacuum energy contribution. We find that the occurrence of phenomenologically interesting minima restrict the 
vacuum energy to lie within a rather restricted range at the EW scale, translating into a corresponding restricted 
range A° ac ~ (3 — 10)m^ 2 , at the GUT scale, when taking into account present LHC and other phenomenological 
constraints. The main practical consequences for phenomenology is to provide additional constraints on top of standard 



34 




FIG. 12. Non thermal relic density in the (mi/2,tan/3) plane for mo = .5mi/j, Ao = 0. Same captions as in Fig. 8. 



mSUGRA parameter constraints, due to the tight connexion between t]q and non trivial miw minima, as well the 
m i/2 — Bo~ tan (3 correlations from Bq input within the no-scale framework. Allowed regions are very restricted when 
considering the strict no-scale boundary conditions with m Q = A Q = 0, providing theoretical exclusion domain prior 
to any additional experimental constraints. 

Concerning the dark matter relic density, a considerably enlarged allowed region of the mSUGRA parameter 
space can be obtained provided one assumes the gravitino to be the true LSP, accounting for the observed relic 
density with important thermal contributions. Perhaps of particular interest is the fact the strict no-scale model 
Bq = mo = Ao = is not excluded by present LHC and other experimental constraints; in particular there is a range 
for toi/2 ~ 400 — 800 GeV, tan/3 ~ 20 — 25 compatible with to/, > 114 GeV, the g^ — 2 deviation (Fig. 7), provided the 
reheating temperature is 10 s — 10 9 GeV as illustrated in Fig. 9. Incidentally this is rather close to the (m^, tan/3) 
range also preferred in recent analysis of the flipped- SU(5) no-scale scenario [25]. A direct comparison of our results 
with the ones in these papers is however limited, since the flipped SU(5) model is quite different, with modified RGEs 
affecting in particular the runnings in the gaugino sector, allowing for a neutralino LSP even for larger TO1/2 than 
found in our study. 

Even if departing slightly from the original no-scale scenarios, the idea of dynamically fixing the soft breaking 
masses from extra minimization at the EW scale remains very attractive, even more so as it emphasizes the role of the 
vacuum energy being crucial for the occurrence of non-trivial no-scale minima. Any future experimental determination 
or exclusion of muj interpreted within no-scale supergravity framework will thus help pinpoint information related 
to the vacuum energy contribution at the EW and possibly at the GUT scale. 



[1] E. Witten, Nucl.Phys. B188, 513 (1981). 
[2] N. Sakai, Z.Phys. Cll, 153 (1981). 



35 




TOi/2 (GeV) 



FIG. 13. Dark matter relic density values consistent with WMAP constraints(5.8) in the (mi/2, tan/3) plane, calculated 
assuming a gravitino G LSP in most of the parameter space, or otherwise indicated, mo = .5m 1 ^ 2 , Ao — and other captions 
as in Fig. 7. The different levels of red corresponds to different reheating temperatures, as indicated. 
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